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ABSTRACT 

We have measured the orbital parameters of seven close binaries, including six new objects, 
in a radial velocity survey of 38 objects comprising a hot subdwarf star with orbital periods 
ranging from ~ 0.17 to 3 d. One new system, GALEX J2205-3141, shows reflection on a M 
dwarf companion. Three other objects show significant short-period variations, but their or¬ 
bital parameters could not be constrained. Two systems comprising a hot subdwarf paired with 
a bright main-sequence/giant companion display short-period photometric variations possibly 
due to irradiation or stellar activity and are also short-period candidates. All except two can¬ 
didates were drawn from a selection of subluminous stars in the Galaxy Evolution Explorer 
ultraviolet sky survey. Our new identifications also include a low-mass subdwarf B star and 
likely progenitor of a low mass white dwarf (GALEX J0805-1058) paired with an unseen, 
possibly substellar, companion. The mass functions of the newly identified binaries imply 
minimum secondary masses ranging from 0.03 to 0.39 Mq. Photometric time series suggest 
that, apart from GALEX J0805-1058 and J2205-3141, the companions are most likely white 
dwarfs. We update the binary population statistics: Close to 40 per cent of hot subdwarfs have 
a companion. Also, we found that the secondary mass distribution shows a low-mass peak 
attributed to late-type dwarfs, and a higher-mass peak and tail distribution attributed to white 
dwarfs and a few spectroscopic composites. Also, we found that the population kinematics 
imply an old age and include a few likely halo population members. 

Key words: binaries: close - binaries: spectroscopic - subdwarfs - white dwarfs - ultraviolet: 
stars. 


1 INTRODUCTION 

Hot subdwarf stars (see a review by |Heber|2009| ) are core helium 
burning stars with very thin hydrogen envelopes and belong to the 
extreme horizontal branch (EHB). The mass of most hot subdwarfs 
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is about 0.5 Mq. The origin of EHB stars, i.e, the hot, hydrogen- 
rich (sdB) and helium-rich subdwarf (sdO) stars, is closely linked 
to binarity. fMengel, Norris & Gross| ( [T976| l first proposed that sdB 
stars are formed in close binary systems and [Dorman, Rood, &| 
|0’Connell| (|1993f inferred the presence of an extremely thin hy¬ 
drogen envelope (< 0.001 Mq). [Han et ah] ( |2002||2003| ) proposed 
three formation channels for sdB stars through binary interaction, 
i.e., common envelope (CE), Roche lobe overflow (RLOE), and bi¬ 
nary merger. [Han et al.| ( |2003) predict a binary fraction of 76 - 89 
per cent with orbital periods ranging from 0.5 hr to 500 d. How¬ 
ever, they caution that the observed frequency could be much lower 
due to selection effects. The proposed formation channels also pre¬ 
dict single sdB stars that form via the merger of two helium white 
dwarfs. Approximately 11-26 per cent of subdwarfs are expected 
to form via this merger channel jHan et al.|2003| >. 

Formation channels of helium-rich (He-sdO) stars are not as 
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well defined. pustham et al.| ( |20lT) proposed that these objects may 
form in a close double degenerate binary with the massive compo¬ 
nent accreting from a helium white dwarf companion and initiating 
helium-shell burning. A small number of sdO stars are known to 
exist as companions to Be stars (|Gies et al.|1998||Peters et al.|2008| 
|2013| l. These sdO stars are formed through close binary interaction 
where the more massive primary star begins mass transfer onto its 
less massive companion during its shell-hydrogen burning phase. 
The result of this mass transfer leaves a spun up Be star with an 
sdO companion ( |Pols et al.|1991| >. 

Cool companions to hot subdwarf stars can be revealed as in¬ 
frared excess in the spectral energy distribution (SED). |The j ll et] 
( |1995^ and |Ulla & Thejll| |T9^ detected infrared excess in 
over 20 per cent of the hot subdwarf stars studied in their sam¬ 
ple. IGirven et al.| ( |2012^ explored photometric surveys that cover 
a wide wavelength range, from the Galaxy Evolution Explorer 
(GALEX) ultraviolet survey through to the infrared, the Two Mi¬ 
cron All Sky Survey (2MASS) and the UKIRT Infrared Deep Sky 
Survey (UKIDSS), and searched for main-sequence companions to 
hot subdwarf stars. They found that the most common companions 
to hot subdwarfs have a spectral type between FO and KO, while 
M-type companions were found to be much rarer. 

Radial velocity surveys (e.g., Maxted et al. 2001 JMr)^es-| 


|Rueda et al.|2003||Copperwheat et al.|2011[ Geier et al. 201 la| ( of 

sdB stars have shown that approximately half of all sdB stars re¬ 
side in close binary systems with either a cool main-sequence star 
or a white dwarf companion. These surveys target binary systems 
with periods of a few hours to 30 days. |Napiwotzki et al.|p004) 
reported a binary fraction of 39 per cent of sdB stars from the 
ESO Supernovae type la Progenitor surveY (SPY). [Copperwheat| 
|et al.|p0TT| > estimated a higher binary fraction of 46 - 56 per cent 
from their survey of sdB stars selected from the Palomar-Green and 
Edinburgh-Cape surveys. 

A few rare sdB stars are found in close orbit with a mas¬ 
sive white dwarf (Mwd <: 0.9 Mo), making them Type la super¬ 
nova progenitors. These systems would first evolve to AM CVn 
systems before detonating either as a Type la or the less ener¬ 
getic Type .la (lax) supernova jBildsten et al.|[2()07[ [Fink et ^ 
|2010 Solheim|2010^ . The first such candidate is KPD 1930-1-2752 
( [Maxted et al.|2000a[ [Geier et al.|2007[ l, with a second candidate, 
GALEX J1411-3053 (CD-30 11223), discovered as part of our ra¬ 
dial velocity survey of GALEX selected hot subdwarf stars (|Vennes| 
[et al.[2012[ (. 


Some sdB stars in close binary systems have stellar param¬ 
eters that fall below the zero-age horizontal branch and probably 
did not initiate helium burning. Such objects have very low masses 
(sj 0.2 Mq) and are the progenitors of extremely low mass (ELM) 
white dwarfs, which will in time evolve into AM CVn systems. If 
the companion to these low mass stars is a massive enough white 
dwarf, then the system may become a Type la supernova. The first 
known low mass sdB star, HD 188112, was discovered by|Heber et| 

[^ j200^ . 

[Ahmad, Jeffery, & Fullerton] ( [2004^ discovered the first dou¬ 
ble subdwarf binary, PG 1544-1-488. This helium rich sdB (He- 
sdB) binary remains, at the present time, unique. The mass ratio 
determined from the velocity semi-amplitude of the components 
show that they have a similar mass which suggests that the sys¬ 
tem emerged from a CE comprised of two nearly identical red gi¬ 
ant cores ( [§ener & Jefrery|2()T4l (. Alternatively, [Lanz et al.|p004) 
interpreted the peculiar atmospheric composition of He-sdB stars, 
such as PG 1544-r488, with evolutionary models involving a de¬ 
layed helium-core flash and convective mixing while descending 


on the white dwarf cooling track. Similarly, HE 0301-3039 is a 
close binary consisting of two sdO stars (jLisker et al.|2004||Stroeer| 
[et al.|2007|l that may be the outcome of double-core CE evolution 
( [Justham et al.|201 1^ . 

Surveys of hot subdwarfs involving photometric time series 
have uncovered several more low mass sdB stars. Kepler obser¬ 
vations revealed that KIC 6614501 is another low mass sdB plus 
white dwarf system ( [Silvotti et al.[2012^ . Also, |Maxted et al.H2014[ l 
presented 17 eclipsing systems from Wide Angle Search for Planets 
(WASP) survey that are likely to contain a pre-helium white dwarf, 
similar to the system ISWASP 1024743.37-251549.2 ([Maxted et[ 
[aTjllTO) . Follow-up spectroscopy for six of these systems con¬ 
firmed them to be main-sequence A stars with very low mass 
(x 0.2 Mq) pre-He white dwarfs currently experiencing hydrogen- 
shell burning. 

Wider binaries (orbital periods ~ years) containing a sdB star 
with a cool main-sequence companion were reported by [Barlo^ 
[et al.[ ( [2012[ [2013a^ and [Vos et al.[ P013[ l. The predicted period 
distribution by Han et al^ i 2003 l is bimodal with some B to F type 
companions in the longer-period range: The relative frequency of 
short- to long-period binaries depends on the actual value of the 
critical mass ratio for stable mass transfer; this ratio may be set 
with a study of potential subdwarf plus A-star binaries. [Chen et[ 
[ar] ( [2013} showed that these long period binaries are the result of 
stable RLOF. 

Venues, Kawka & Nemeth| j2011[l and [Nemeth, Kawka &| 


[Vennes[ ( [2012[ l presented a new sample of sdB stars selected 
from the GALEX all-sky survey and we conducted a radial ve¬ 
locity survey of a subsample of stars from this selection. The 
first two systems (GALEX J0321-H4727 and GALEX J2349-(-3844) 
discovered as part of this survey were presented by ||Kawka et| 
[aTI pOTO) , followed by the aforementioned short-period system 
GALEX 11411-3053 ( [Venues et ar]|20T2l [GeTer et al.[[2013a[ l. 
Additional spectroscopic and photometric observations of the first 
two systems were presented in [Kawka et al.[ ( [2012a[ l along with 
a progress report on the other systems that were observed as part 
of this programme. The photometric observations confirmed the 
reflection effect in GALEX 10321 -h 4727 originally reported by 
[Kawka et al.[j20ro) and based on Northern Sky Variability Sur¬ 
vey (NSVS) photometry. The observations also showed that both 


GALEX 10321-t4727 and GALEX 12349-(-3844 are V2093 Her 
type pulsating subdwarfs ( [Green et al.|2003^ . 

In this paper, we present spectroscopic and photometric ob¬ 
servations of a sample of GALEX-selected hot subdwarf stars with 
the aim of determining their binary properties. Sections 2.1 and 2.2 
present details of our spectroscopic observations, while Section 2.3 
present archival photometric time series. In Section 3 we present 
an analysis of stellar properties (3.1), and of binary properties sup¬ 
plemented by our analysis of photometric time series (3.2). Einally, 
we present a review of the properties of known binaries compris¬ 
ing a hot subdwarf star, including the properties of the components 
(Section 4.1), the population kinematics (4.2), and the properties of 
some outstanding individual cases (4.3), followed by a summary of 
the present work (4.4). 


2 SAMPLE SELECTION AND OBSERVATIONS 

Table[T]lists the stars originally included in our radial velocity sur¬ 
vey with notable properties described in Section 3.1. The sample 
includes 38 spectroscopically confirmed hot subdwarf stars, and 
two objects that were respectively identified as an early B star 


© 2015 RAS, MNRAS 000,[Tp6[ 

















































































































Hot subdwarf binaries 3 


Table 1. Target summary. 


GALEX J 


Other names 


Tea logg log(He/H) Notes" 

(K) c.g.s. 


004759.6+033742 
004729.4+095855 
004917.2+205640 
005956.7+154419 
020656.1 + 143900 
023251.9+441126 
040105.3-322348 
050018.9+091203 
050735.7+034814 
061325.3+342053 
065736.7-732447 
070331.5+623626 
071646.9+231930 
075147.1+092526 
080510.9-105834 
081233.6+160123 

104148.6- 073034 
111422.0-242130 

135629.2- 493403 
140747.6+310318 
141133.3+703737 
142126.5+712427 

142747.2- 270108 
143519.8+001352 
163201.4+075940 
173153.7+064706 
173651.2+280635 

175340.5- 500741 
184559.8-413826 

190211.7- 513005 
190302.4-352828 

191109.2- 140651 

203850.3- 265750 

215340.4- 700430 

220551.8- 314105 
225444.1-551505 
233451.7+534701 

234421.6- 342655 


BPS BS 17579-0012, PB 6168 
HD 4539, HIP 3701 
PG 0046+207 

HIP 4666, PG 0057+155, PHL 932 

CHSS 3497 

FBS 0229+439 

CD-32 1567, EC 03591-3232 

HS 0457+0907 


CPD-73 420 
FBS 0658+627 
TYC 1909-865-1 

TYC 5417-2552-1 

TYC 5492-642-1 

EC 11119-2405, TYC 6649-111-1 
CD-48 8608, TYC 8271-627-1 
BPS BS 16082-0122 
TYC 4406-666-1 
TYC 4406-285-1 

EC 14248-2647, TYC 6740-942-1 
TYC 325-452-1, PG 1432+004 
TYC 960-1373-1, PG 1629+081 

TYC 2084-448-1 


CD-51 11879, TYC 8386-1370-1, LSE 263 
BPS CS 22936-0293 
TYC 5720-292-1 
TYC 6916-251-1 

EC 21494-7018, TYC 9327-1311-1 
TYC 7489-686-1, BPS CS 30337-0074 

TYC4000-216-1 

CD-35 15910, HE 2341-3443 


33260^1“ 

30490+^*” 

-23,0 

36270+f®Q 

2-i99Wf° 

34250;|« 

29940^® 

28750!3™ 

11140/9310 

30620/^0 

22320+^30 

31580/™ 

27440/|“ 

23430/™ 

33070+^"® 

24900+^° 

21170/;^“ 

25620/™ 

31880/™ 

23090/ 

38110/™ 

27780/'™ 

36160/|“ 

32430+sSo 

-^70 

72300 

32100+^"^o 

55970+^5^0 

jj^/u_. 7 gQ 

58450/™“ 

23720/™ 

28650/™ 

31070/™ 

35680/™ 

28390/'*'“ 


6 14+0-22 
c oo-f0.03 

c OO+0.04 

-Q-P5 

c 77+0.05 

-Q-Q6 

c 7'3+O.Oy 
c 71 +0.06 

75^8:?^ 

''^•''^-0.13 

-r 47+0.08 

'^-Q-Q6 
^ 4 C+ 0.07 

4o78:45 

-0.04 

4.39/3.67 

5 74+0 -11 
- 0.12 

c ^o+0.03 
^•°°-0.06 
5.56/ 
5.63/ 

5 29+^-^^ 
'^•"^^-0.07 
-r 74+0.07 
-0.16 
4 7 C+O .03 
^•^'^-0.09 


m 


0.23 

0.04 


5.67: 

5 70+0-05 
-Q-Q8 

c 70+O.O8 
^•^^-0.08 
c 00+O.O6 
-0.09 

<r o<r+0.I8 

^•■^•^-0.07 
^ 74+0.84 

59^8:ft 

-Q-IS 
c 772+0.27 

26^8^1 

5 69+0-^i 
-Q-09 
5 04+° 

^•^'^-0.02 
c ^o+O.OT 
^•^°-Q.Q3 

^•^^-0.06 
5 01 +0.07 
^•^^-0.06 
c oq+0.05 


:8:4i 


-2 63+044 

_7 47 + 0.20 
-0.07 
-2 48+046 

-169^8:8^ 
_ +8:?! 

^■°'-0.24 
-1.70 

-1.92 
-146™ 
-305-8:i^ 

-Q-7§ 
1 70+0.04 

^-40-0 08 

< -3.21 

-2 76+0-22 
-0,26 

_7 49+0.27 
^•^^-0.30 

< -3.44 

< -2.90 

-9 44+0 16 

-246-8:?? 

-0.31 
—9 7S+0-25 
^•^'^-0.43 
-1.18+““* 


-0.09 

-2.36 

-3.7 


< ■ 

< 

-1.71/ 

-2.39+ 


-2 71+ 

-0.29 

< -2.53 

0+0.69 


8;?9 


hit 


-1.09; 

-2.25 

+2 10+*:'''® 
-0-38 

+0 02 +”'“ 
+U.UZ jjuj 

< -1.96 
+0 25+“ ’'“ 

, n;8:§9 

'■*■^-029 

77+0.13 

-’■^^-1.15 

-2 09+“ 
-Q93 
-2 47+0 '= 
-0.13 

-1.43 ±0.07 

-3 07+“'^' 
-’+''-0.26 


sdB+F6V; IR; nearby star 


CaH&K, RV 
RV 

nearby star 

close B+A V binary, RV 
nearby star (6 arcsec), RV 
ELM WD progenitor, RV 
RV 


sdB+G8V; IR 
high-o early B 

sdB±F; IR; ELM WD progenitor? 


nearby star, RV 
RV 

sdB+F7V; IR; variable 
sdB+F7V; IR 
sdO; He I spectrum 
sdO; Hen spectrum 
RV 

sdO; Hen spectrum 
sdO±G3.5III; IR; variable 
ELM WD progenitor? 
reflection, RV 
RV 



Other names 

Teff 

logg 

log(He/H) 

Notes 

I 


(K) 

c.g.s. 



123723.5+250400 

Feige 66 

34300+™ 

34640/™ 

5.82 ± 0.04 

-1 51+0 05 

l.Dl 


160011.8-643330 

TYC 9044-1653-1 

6 02+0 08 
O.u/,_o 11 

-0 30+0 05 
-0.04 



RV: confirmed radial velocity vaiiable star; IR: SED of the stars shows significant IR excess. 


and a AV+BV binary. The early B star GALEX JI407+3103 
is notable for its high radial velocity, while the close AV+BV 
pair GALEX J0716+2319 shows significant radial velocity vari¬ 
ations on a short time scale. All except two objects were ran¬ 
domly selected from our catalogue of GALEX/Guide Star Cata¬ 
logue ultraviolet-excess obje cts jVennes, Kawka & Nemeth|201H 
[Nemeth, Kawka & Vennes||2012| . Briefly, the source catalogue 

includes bright objects (Nuv < 14) with an ultraviolet excess 
(Nuv - V < 0.5). The latter criterion still allows for the selec¬ 
tion of hot subdwarf plus F/G dwarf pairs (see [Vennes, Kawka| 
|& NemelhllTOllI ). Two additional stars that were not observed 
by GALEX, including a blue-excess object jJimenez-Esteban, Ca-| 
[ballero, & Solano|201 1^ , are listed at the bottom of Table [T] with 
J2000 coordinates. 


The GALEX name corresponds to the coordinates of the ul¬ 
traviolet source detected in the near ultraviolet (NUV) band (Sec¬ 
tion 2.3); for convenience, the names are abbreviated to 4 dig¬ 
its right ascension and declination. The ultraviolet coordinates 
are generally close to the Guide Star Catalog (GSC2.3.2) opti¬ 
cal coordinates (<1 arcsec), but, in a few cases, offsets as large 
as 4 to 9 arcsec occurred (GALEX J1421+7124, J1427-2701, 
J1902-5130, J2344-3426). Despite the offsets, the ultraviolet and 
optical sources must be one and the same. These offsets cannot be 
attributed to a high proper-motion and are most likely due to a dis¬ 
torted point spread function (PSF) in bright off-centred sources in 
the GALEX images (Section 2.3). 

Throughout this paper we will refer to the hot subdwarf as the 
primary and its companion as the secondary. Table[T]lists some no¬ 
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table particularities such as the presence of a nearby star, whether 
unrelated or physically associated to the hot subdwarf, a bright 
main-sequence companion, or photometric variability due to reflec¬ 
tion on a late-type companion or stellar activity (see Section 3.1). 
Most stars display Hi-dominated line spectra, but we also noted 
the presence of He-rich subdwarfs characterized by He i and He ii- 
dominated line spectra. The stellar parameters of a handful of sub¬ 
dwarfs locate them below the zero-age EHB (ZAEHB) and these 
objects are likely progenitors of ELM white dwarfs (Sections 3.1 
and 3.2). 

2.1 Intermediate to high-dispersion spectroscopy for radial 
velocity measurements 

Our first extensive set of observations was obtained with the Wide 
Field Spectrograph (WiFeS, [Dopita et al.||2007} attached to the 
2.3 m telescope at the Siding Spring Observatory (SSO). The ob¬ 
servations were conducted on UT 2011 July 14 to 18, UT 2011 De¬ 
cember 2 to 3 and UT 2012 April 27 to 30. We used the B3000 and 
R7000 gratings with a slit width of 1 arcsec that provided spec¬ 
tral ranges of 3200-5900 A at a resolution of R = T/A/i = 3000 
and 5300-7000 A at R - 7000, respectively. The RT560 dichroic 
beam splitter separated the incoming light into its red and blue com¬ 
ponents. WiFeS is an image-slicing spectrograph with 25 slitlets 
(38 X 1 arcsec) and depending on the seeing, the target can cover 
a few slitlets. The signal-to-noise ratio (S/N) of each observation 
was maximized by extracting the spectrum from the most signif¬ 
icant (< 6) traces. Each trace was wavelength and flux calibrated 
prior to co-addition. The spectra were wavelength calibrated using 
NeAr arc spectra that were obtained either prior to or following 
each observation. 

Next, our second set of observations was obtained using the 
Ritchey-Chretien Focus (R.-C.) Spectrograph attached to the 4 m 
telescope at Kitt Peak National Observatory (KPNO) on UT 4 - 6 
January 2012. We used the KPC24 grating in second order com¬ 
bined with the T2KA CCD to provide a spectral range of 6030 - 
6720 A and a dispersion of 0.52 A pixeU*. The slit width was set 
to 1.5 arcsec which provided a resolution of ~0.9 A or R = 7000. 
Contamination from third order was removed using the GG495 fil¬ 
ter. The spectra were wavelength calibrated using HeNeAr spectra 
which were obtained following each observation. 

We obtained a third set of observations using the ESO Faint 
Object Spectrograph and Camera (EFOSC2) attached to the 3.6 m 
New Technology Telescope (NTT) at La Silla Observatory in 
September 2012. We used grism number 20 centred on Ho- pro¬ 
viding a spectral range from 6040 to 7140 A and a dispersion of 
0.55 A pixeL'. We set the slit width to 0.7 arcsec resulting in a 2 A 
resolution or R = 3500. Next, we obtained additional spectra with 
EFOSC2 on the NTT on UT 31 July and 1 August 2014. We used 
grism number 19 that provided a spectral range from 4435 to 5120 
A and, after binning 2 x 2, a dispersion of 0.67 A per binned pixel. 
The slit width was set to 1 arcsec resulting in a resolution of ~2 A 
or R ^ 2000. Additional EFOSC2 spectra of GALEX J1731+0647 
were extracted from the ESO archive (programme 090.D-0012, PI 
S. Geier). The data were also obtained with grism 19, but binned 
2x1 resulting in a dispersion of 0.34 A pixeL*. The slit width was 
set to 1 arcsec resulting in a resolution of ~2 A. All spectra were 
wavelength calibrated using HeAr arc spectra which were obtained 
following each observation. 

Also, we obtained a fourth set of spectra using the grating 
spectrograph attached to the 1.9 m telescope at the South African 
Astronomical Observatory (SAAO) on UT 2014 February 11. We 


used the 1200 lines/mm grating with a blaze wavelength of 6800 
A. This arrangement provided a range of 6023 to 6782 A with a 
dispersion of 0.439 A per pixel. The slit width was set to 1.05 arc¬ 
sec resulting in R = 7000, or a resolution of «1A at Hu. A CuNe 
comparison arc was obtained following each target observation. 

We assembled a fifth data set with observations of the bright 
objects HD 4539 (GALEX J00474-0958), the spectro-photometric 
standard Feige 66, GALEX J1421-)-7124, GALEX J1736-I-2806, 
and GALEX J23344-5347 using the 2m telescope at Ondfejov Ob¬ 
servatory. The observing configuration and procedure are described 
in |Kawka et al.| ( |2010^ . Briefly, for each star we obtained a series 
of spectra centred on Ha. We used the 830.77 lines mm“* grating 
with a SITe 2030 x 800 CCD, this resulted in a spectral resolution 
of R = 13 000. Each target exposure was immediately followed by 
a ThAr comparison arc. 

Finally, and introducing our sixth and most recent observation 
programme, we obtained three high-dispersion echelle spectra of 
the short period binary GALEX J2254-5515. From UT 24 Novem¬ 
ber to 4 December 2014 we used the Fiber-fed Extended Range 
Optical Spectrograph (FEROS) attached to the 2.2 m telescope at 
La Silla. The spectra range from =3600 to »:9200 A at a resolution 
ofR^ 48,000. 

We supplemented our data sets with archival spectra. We ex¬ 
tracted processed FEROS data from the ESO archive. The spec¬ 
tra were obtained under the programmes 076.D-0355, 077.D-0515, 
078.D-0098 (PI: L. Morales-Rueda) and 086.D-0714 (PI: S. Geier). 

We also extracted spectra from the Isaac Newton Group 
(ING) Archive. The first set of data (GALEX J1632-H0759 and 
GALEX J1731-)-0647) was obtained with the Intermediate Disper¬ 
sion Spectrograph (IDS) attached to the Isaac Newton Telescope 
(INT) on UT 17 May 2013 (run numbers 984456 and 984458) and 
on UT 19 May 2013 (run numbers 984760, 984762 and 984763). 
The spectra were obtained with the R1200B grating which resulted 
in a useful range of 3900 to 5200 A and a dispersion of 0.48 
A and delivering a resolution of 1.5 A assuming a 3-pixel full- 
width at half-maximum (FWHM). The spectra were wavelength 
calibrated using CuAr and CuNe arcs and adjacent exposures were 
co-added to obtain the final radial velocity. A second set of data 
(GALEX J1632-I-0759) was obtained with the William Herschel 
Telescope (WHT) and the Intermediate dispersion Spectrograph 
and Imaging System (ISIS) on UT 26 August 2010 (run numbers 
1483813 and 1483814). The spectra were obtained with the R600B 
and R600R gratings and calibrated with CuAr and CuNe arcs re¬ 
sulting in useful ranges of 3500-5100 A and 5500-7030 A and dis¬ 
persion of 0.88 A per binned pixel in the blue (2 x 2) and 0.49 A 
pixeL* in the red (binned 2 x 1), corresponding to spectral resolu¬ 
tions of 1.7 A in the blue and 1.5 A in the red assuming a 3-pixel 
FWHM. 

On average, a high signal-to-noise ratio was achieved with 
EFOSC2 on the NTT (^ ^ 100), the R.-C. Spectrograph on 
the KPNO 4 m telescope (^ ^ 60), and WiFeS on the SSO 
2.3 m telescope (S/N » 80). A lower signal-to-noise ratio was 
achieved with the coude spectrograph on the Ondfejov 2 m tele¬ 
scope, FEROS on the MPG 2.2 m telescope (La Silla), and the 
grating spectrograph on the SAAO 1.9 m telescope (S/N 30). 
The lower S/N achieved at Ondfejov and La Silla is largely com¬ 
pensated by the higher dispersion resulting in comparable or supe¬ 
rior velocity accuracy (see next Section). More than 70 per cent of 
our spectra had a S/N> 40 and spectra with ill-defined hydrogen or 
helium lines (S/N< 15) were rejected. 
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Feige 66 




TYC 4000-216-1 




Figure 1. Line profile analysis of the WHT spectrum of Feige 66 (top) and KPNO spectrum of TYC4000-216-1 (bottom), labelled with best-fitting parameters. 


2.1.1 Tests of the wavelength and velocity scales 

We performed a series of tests of the wavelength scale of relevant 
spectra using the O i sky emission lines and atmospheric molecular 
absorption bands. 

Diffuse 01/16300.304 emission helps set the accuracy of the 
wavelength scale, particularly in low- to intermediate-dispersion 
spectra. A strong emission line is detected in 93 per cent of 
all usable EFOSC2 spectra, 95 per cent of all KPNO and SSO 
spectra, and nearly all SAAO spectra. A short exposure time as 
well as the appearance of scattered moonlight usually limit the 
usefulness of this template. The Oi velocity averaged d(Oi)= 
O.Okms^* at KPNO with a dispersion (T„(Oi)= 2.1 kms^*, d(Oi)= 
1.9kms“* with EFOSC2 and a dispersion cr„(Oi)= 5.4kms“*, 
i;(0 1 )= 3.7 km s“' at SSO and a dispersion cr„(0 1 )= 7.4 km s“*, and 
i;(Oi)= 4.6kms“* at SAAO and a dispersion cr„(Oi)= 4.4kms“*. 


The emission line appeared blended in most spectra obtained dur¬ 
ing bright time at SSO. Based on this analysis the expected ac¬ 
curacy should be of the order of 2-7 km s“'. The accuracy of the 
wavelength scale using coude or echelle spectrographs is normally 
of the order of 1 km s“' or better. 

Systematic velocity shifts are expected following an improper 
placement of the star on the slit, particularly if the stellar image is 
much narrower than the slit width. Excellent seeing conditions are 
often encountered at La Silla and KPNO. We cross-correlated tel¬ 
luric absorption features in the KPNO and EFOSC2 spectra with a 
telluric template of identical spectral resolution. We measured an 
average velocity of-0.8kms“' with a dispersion of 13.4kms“' in 
the EFOSC2 spectra and an average velocity of 2.4kms“* with a 
dispersion of 10.0 km s“* in the KPNO spectra. Velocity deviations 
of up to 50 km s“* were found in a few well-exposed EFOSC2 spec- 
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tra: We corrected the measured stellar velocities at La Silla using 
the telluric template velocities. 

In summary, after applying telluric corrections, we estimate 
that errors in stellar velocity measurements due to various sys¬ 
tematic effects are better than ~10kms“' provided that the pho- 
tospheric lines are well defined. Ultimately, the accuracy of the 
wavelength is verifiable using actual stellar data and by plotting 
the velocity dispersion distribution (Section 3.1.3). 


2.2 Low-dispersion spectroscopy for stellar parameter 
determinations 

For stars not listed in [Nemeth, Kawka & Vennes|P012] l, we ob¬ 
tained additional low dispersion spectra with the R.-C. spectro¬ 
graph attached to the 4m telescope at KPNO on UT 2013 July 12 
(GALEX J1421+7124) and 2014 May 24 (GALEX J2334-)-5347). 
We used the KPC-lOA grating and T2KA CCD with a dispersion 
of 2.77 A pixeL’ in first order and centred on 5875 A. We used the 
order sorting filter WG360 and set the slit width at 1.5 arcsec result¬ 
ing in a spectral resolution of x: 5.5 A. The spectra were wavelength 
calibrated using the HeNeAr arc. 

We extracted a set of spectra of the spectro-photometric stan¬ 
dard Eeige 66 from the ING archive. These spectra were obtained 
with ISIS attached to the WHT (run numbers 133198, 133200, 
133201, 133223, 133226, 133227). The spectra were obtained us¬ 
ing the R300B grating in the blue arm providing a dispersion of 
1.54 A pixeL' and a spectral range from 3620 to 5190 A. The slit 
width was set to 2.4 arcsec for each observation which corresponds 
to a resolution of x: 7.5 A. The spectra were wavelength calibrated 
using a CuAr arc. 

Details of the low-dispersion spectroscopy obtained of other 
objects in the present sample are given by Vennes, Kawka &| 
|Nemeth|j2011[ l and |Nemeth, Kawka & Venn^ \2012) . 


2.3 Photometry and imaging 


We compiled available optical and infrared photometric measure¬ 
ments and combined them with the GALEX NUV and EUV pho¬ 
tometry from the all-sky imaging survey (AIS) to build a SED for 
each object in the sample. The ultraviolet data were collected from 
the site galex. stsci .edu/GalexView/. [Morrissey et'dl|j2007f 


present details of the instrument calibration. Table Al in the Ap 


pendix lists the GALEX magnitudes, along with the available V 
magnitudes as well as 2MASS ( [Skrutskie et al.|2006 i and Wide- 
field Infrared Survey Explorer (WIS E, Wright et al.|2010^ infrared 
measurements. 

The PSF of WISE images ranges from 6 to 12 arcsec in the 
3 to wavelength range, while the PSF in 2MASS images is 
close to 2.5 arcsec. Because of its relatively broad PSF, stars lo¬ 
cated within its range and identified in higher-resolution imaging 
are certainly contaminating the SED in the mid-IR range. 

Also, we extracted photometric time series from the Su- 
perWASP (SWASP, [Mlacco et al.||2006) public archive, NSVS 


I Wozniak et al.|2004 


_ , All Sky Automated Survey (ASAS;|Poj man-| 

|ski||1997^ and Catalina surveys ( [Drake et al.||2009^ . The Catalina 
photometry is unfiltered. Bright targets (< 12 mag) are often sat¬ 
urated, but the photometric measurements are more precise with 
faint targets (> 14 mag) than those obtained in the other three sur¬ 
veys consulted. The SWASP images are filtered (4000-7000 A). 
Light curve analysis of SWASP data is valuable because of the 
large number of measurements obtained for individual targets. We 


obtained AS AS times series in the V band and the NSVS images 
are unfiltered. 

The calibrated GALEX magnitudes are obtained from the 
count rates extracted using elliptical apertures (fuv_flux_auto, 
nuv_flux_auto) fitted to the actual stellar profiles and converted 
into the AB system. The average GALEX PSP is matched approx¬ 
imately by Gaussian functions with FWHM of 5.3 and 4.2 arc¬ 
sec in NUV and FUV images, respectively, and a positional ac¬ 
curacy of X! 0.5 arcsec. However, several factors affect the relia¬ 
bility of the GALEX photometric magnitudes. The GALEX imag¬ 
ing quality varies with the detector position with a strong depen¬ 
dency on the radial distance from the image centre. We recorded 
the target distance to the centre of the field of view (fov.radius), 
as well as the actual FWHM values in the FUV and NUV im¬ 
ages (fuv j'whm.world, nuvTwhm.world) for each target. Mea¬ 
surements with a radial distance outside of 0°.4 combined with a 
large PSF (> 0°.01) or measurements with an exceedingly large 
PSF (> 0°.04) are marked in Table [aT] as possibly unreliable. Fi¬ 
nally, bright objects with unreliable non-linearity corrections out¬ 
side the range of validity are marked. Non-linearity effects dom¬ 
inate the photometric error: A 10 per cent loss is observed at 
N^v = 13.9 and F^v = 13.7 so that most measurements in the 
present selection are affected. [Morrissey et al.[([2007| and[Camarota[ 
[& Holber^ 1 2014) propose correction algorithms that are nearly 


Camarota & Holberg[([2014]( presented a calibration s 


identical. 

pie sufficiently large to allow us to evaluate the scatter in the syn¬ 
thetic versus measured magnitude relations. For example, this scat¬ 
ter is of the order of 0.35 and 0.4 mag at Nyy and Euv = 13, 
respectively. We adopted GALEX magnitudes adjusted using the 
correction algorithm of [Camarota & Holberg| ( [2OT^ with errors es¬ 
timated using the scatter in these corrections for a given magnitude. 

The EFOSC2 acquisition images provide a deep, high-spatial 
resolution view of the fields surrounding target stars. These images 
were obtained with the Loral/Lesser 2048x2048 CCD. With a focal 
plane scale of 8.6 arcsec mm“' and a pixel size of 15 /rm, the sky 
images are sampled with a pixel size of 0.129 x 0.129 arcsec^, or, 
after binning 2x2, a pixel size of 0.258x0.258 arcsec^. The images 
allow for the identification of physical companions or unrelated, 
nearby stars. 

Fig. [Al[ Fig. [A2[ and Fig. [A3[ in the Appendix compare all 
available photometry to synthetic spectra computed using stellar 
parameters listed in Table[^ 


3 ANALYSIS 

We present, in order, the properties of the sample including an 
overview of the stellar parameters (Teff, logg, logHe/H) and evo¬ 
lutionary history, the characteristics of the SEDs and photometric 
time series, and the radial velocity data set. We identify new binary 
candidates and present an analysis of individual binary properties 
from the combined data sets. 


3.1 Sample properties 

Table lists the atmospheric parameters obtained from [Nemeth, [ 
[Kawka & Venn^ \2Q\2) . The Balmer line analysis for 
three additional objects (Feige 66, GALEX 11421-1-7124 and 
GALEX J2334-I-5347) is based on the model grids of [Vennes,[ 
[Kawka & Nemeth| pOl l[ l. Best fitting parameters (Teff, logg, 
logHe/H) are obtained usingminimization techniques with the 
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Figure 2. Physical properties, luminosity versus eifective temperature, of 
the sample: sdO stars are shown with full triangles while sdB stars are 
shown with open squares (assuming 0.47 M©) or full squares (assuming 
0.234 Mq). The zero-age EHB is labelled “a” while the terminal-age EHB is 
labelled “b”. Evolutionary tracks computed by |Dorman, Rood, & 0’Connell| 
jl993) with a helium core mass of 0.469M© and hydrogen envelopes of, 
left to right, 0.002, 0.004, 0.006, and 0.01 M© are shown with full lines. 
The cooling track from |Driebe et al.]jl998^ for progenitors of ELM white 
dwarfs of 0.234 M© is shown prior to hydrogen shell flashes with a dashed 
line. Lines of constant radii at 0.01, 0.1, 1, and lOR© are labelled accord¬ 
ingly. 

observed line profiles (Hei,ii and Hi) being simultaneously ad¬ 
justed to interpolated spectra from the model grid. Examples of 
Balmer and helium line analyses are shown in Fig.^ 

Fig. 0 shows properties of the sample presently investigated. 
Using the effective temperature ) and surface gravity (g) we 
determined the total luminosity (in Lq) by adopting for most objects 
a sample-average mass of 0.47 Mg ( [Fontaine et al.|20l^ . 

L = 4nR^a-T^ff, 

where cr is the Stefan-Boltzmann constant and the radius (R) is 
calculated using 



where M is the subdwarf mass and G is the gravitational constant. 

The sdB stars form a sequence of approximately constant 
luminosity, L =10-30 Lgor My = 4.3 (cr = 0.9) mag, and lo¬ 
cated between the ZAEHB and the TAEHB while a few ageing 
sdB stars and all He-rich sdO stars set out on a higher luminos¬ 
ity excursion beyond the stable He-burning stage. The objects ly¬ 
ing below the ZAEHB (L < 10 Lq) with a low-temperature and 
a high-gravity, GALEX J0805-1058 and, tentatively, J1411+7037 
and J2153-7004, were singled-out and were attributed a mass of 
0.23 Mg based on their likely evolutionary status (jDriebe et al.j 
|1998| (. Most objects lie to the left of the EHB tracks suggesting 
that their hydrogen layer is thinner than 0.002 Mg, or, possibly, 
that their surface gravity is overestimated. To investigate the latter 


Figure 3. NUV - V versus V — J colour-colour diagram. Stars with a 
composite IR excess are shown with full black circles and stars with an 
IR excess due to a nearby star are shown with open circles (see Section 
3.1.1), while all others are shown in grey. Models at 24, 28, 32, 36, and 
40 X 10^ K (logg = 5.7, logHe/H = -2.5) are shown, from bottom to 
top, with open squares linked by a full line. The effect of interstellar extinc¬ 
tion (Eb-v = 0.05) on the colours is shown with an arrow in the upper left 
corner. 

possibility, we compared the results of a model atmosphere anal¬ 
ysis using the hydrogen Stark broadening t ables of|Lemke|(|1997|(, 
employed in the present work, to those of jTremblay & Bergeronj 
( [2009^ , which include improved treatment of merging atomic en¬ 
ergy levels. We found that improvements in Stark broadening the¬ 
ory may account for a shift of A log g = -1-0.08 dex near 30 000 K 
(see, e.g., |0stensen et al.|2014[|Telting et al.|2014) in agreement 
with a shift of Alogg = -1-0.06 dex found at 40000 K byjKlepp &| 
|Rauch| ( [2011^ . These systematic shifts are notable, but still cannot 
explain the model offsets apparent in Eig. Metallicity has little 
effect on temperature and gravity below Tefr^ 35 000 K as demon¬ 
strated by jLatour et al.|p014| (. It is worth noting that, while a sdB 
mass of 0.47 Mg may be typical, it does not necessarily apply to all 
objects (e.g., ELM progenitors). 

On the other hand, |Schindler, Green, & Amett| ( [^14[ ( pointed 
out that current evolutionary models fail to reproduce some ob¬ 
served properties of EHB stars, such as the core mass derived from 
asteroseismology, and concluded that evolutionary models must be 
updated to match observed seismic and spectroscopic stellar pa¬ 
rameters. [ScErndjerTGreenr^^m^ ([20^ found that very high 
convective overshooting would be needed to reproduce the seismic 
core mass but that it would, quite improbably, double the EHB life¬ 
time. Therefore, they conclude that the general treatment of convec¬ 
tion in evolutionary models needs updating, and that new opacity 
tables and diffusion calculations are required. 


3.1.1 Overview of the SEDs 

Fig .|^shows the NUV-V versus V- J colour diagram for the sam¬ 
ple of hot subdwarfs listed in TablejAl] Fig.|^and Fig.|^present the 
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J-H 


Figure 4. Same as Fig. I^but showing V — J versus J — H. 



-0.4 -0.2 0 0.2 0.4 

H-Wl 


Figure 5. Same as Fig.|^but showing J-H versus H -Wl. 


14-J versus 7-// and 7-Ff versus H-Wl diagrams, respectively. 
The effect of interstellar extinction is evident in the NUV- V colour 
of some objects, but many are also affected by large systematic er¬ 
rors in the GALEX NUV photometry (non-linearity). For example, 
a larger extinction than observed in the interstellar line-of-sight is 
apparent toward GALEX 116324-0759. The International Ultravio¬ 
let Explorer iWE) spectra that supplement uncertain GALEX pho¬ 
tometric measurements indicate Eb-v = 0.4, largely in excess of 
that found in the extinction map of [Schlegel, Einkbeiner, & Davis] 
( |1998t , Eb-v — 0.08. The effect on ultraviolet colours of an extinc¬ 
tion coefficient Eb-v = 0.05 (see Fig.|^ is relatively modest, but 


coefficients in excess of 0.4 would displace colours from the upper 
left to the reddened, lower right corner in the vicinity of composite 
stars. 


Individual SEDs may be contaminated by the presence of a 
nearby star, either physically associated or unrelated to the hot 
subdwarf. Inspection of the P82, P83, P85 and P89 EFOSC2 ac¬ 
quisition images obtained by Venn es, Kawka & Nemeth] j20I I| l 
and [Ndmefh, Kawka & Venn^ ' j20l2| revealed the presence of 
nearby companions (< 3 arcsec) to 10047-1-0337, J0657-7324, and 
116324-0759 (see below). Visual inspections of the guiding images 
displayed at KPNO did not reveal the presence of a nearby compan¬ 
ion to any other objects. An inspection of photographic plate mate¬ 
rial I http: //surveys. roe. ac. uk/ssa/index. html i helps lo¬ 
cate other, more distant objects (> 3 arcsec) that may contaminate 
photometric measurements with large PSF (e.g., WISE). For ex¬ 
ample, GALEX 107514-0925 is accompanied by a faint (Am ~5 
mag), nearby star ~6 arcsec away at a position angle of 200° (Epoch 
1993). 

The composite nature of the spectra of GALEX 100474-0337, 


114114-7037, J1736-b2806, J1753-5007, 12038-2657 ( |Nemeth,| 
|Kawka & Venne^|2012| ) are confirmed by their IR photometric 
colours. The flux contributions in the V band for these hot subd¬ 
warfs are offset by ~0.7 (possibly contaminated), ~0.0 (weak sec¬ 
ondary detection in the optical), ~L2, ~L3, and ~L3 mag, respec¬ 
tively, relative to their observed composite V magnitudes. Although 
evident in the IR colours of GALEX 11356- 4934 (Eig. | A2^, the 
presence of a companion was not detected by [Nemeth, Kawka ^ 
|Vennes| ( |2012t . The flux contribution from the hot subdwarf in the 

V band is offset by ~0.4 mag relative to the observed composite 

V magnitude. We re-examined the spectra of GALEX 11356-4934 
and we found weak signatures of a cool main-sequence companion 
in the blue spectrum used by [Nemeth, Kawka & Venne^ ( |2012| l, 
and stronger spectral lines representative of a cool main-sequence 
star in the red spectra used in this paper. We present our spectral 
decomposition of this system in Section 3.2.2. The IR colours for 
GALEX 100474-0337, 10657-7324, and 116324-0759 are almost 
certainly contaminated by their nearby, resolved companions. 


3.1.2 Overview of the photometric time series 

Table [^ summarizes the photometric time series analyses. We 
included objects showing significant radial velocity variations 
(e.g., GALEX 12205-3141), objects with composite optical 
spectra (e.g., GALEX 117364-2806), and, finally, three objects 
with previously published analyses from the present survey: 
GALEX 10321-H4727 and GALEX 12349-H3844 ( [Kawka et aLj 
|20T0j (, and GALEX 11411-3053 ( [Vennes et aL]|20T2l (. Photo¬ 
metric variations observed in the sdB plus white dwarf system 
GALEX 11411-3053 are an example of ellipsoidal variations in 
this class of objects. Eourier transform calculations of available 
light curves (Section 2.3) uncovered three objects with signif¬ 
icant periodic variations. The light curves were analysed using 
fast Fourier transform analysis from [Press et ^ ( [1992[ l. Both 
GALEX 117364-2806 and GALEX 12038-2657 are binaries com¬ 
prising a hot subdwarf with a more luminous optical compan¬ 
ion, while GALEX 12205-3141 is composed of a hot subdwarf 
and late-type companion. The photometric variations in the lat¬ 
ter are clearly timed with the orbital period (see Section 3.2) and 
caused by reflection of the primary on the cool secondary. Vari¬ 
ations in GALEX 11736-H2806 and GALEX 12038-2657 may be 
caused by a spot on the surface of the secondary coupled to the 
rotation period. The variable, double-peaked Ha line profile of 
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Table 2. Photometric time series. 


Name 

Survey 

HJD range 
(2450000+) 

Number 

Period range 

(d) 

Semi-amplitude 

(mmag) 

Average magnitude 
(mag) 

Standard deviation 
(mmag) 

100474-0337 

ASAS 

1868-5168 

378 

>0.02 

4.4 ± 7.0 

12.336 

94.5 


NSVS 

1382-1549 

177 

>0.01 

17.8 ± 3.4 

12.676 

32.2 

J0321-I-4727 

NSVS 

1373-1630 

173 

0.26586 “ 

61.3 ±3.9 

12.034 

56.9 


SWASP 

3196-4458 

4575 

0.26586 “ 

43.5 ± 1.0 

11.490 

49.4 

J0401-3223 

SWASP 

3964-4485 

14208 

>0.01 
1.85735 * 

8.4 ±0.1 

0.8 ±0.1 

11.268 

12.4 

J0507-I-0348 

Catalina 

3643-6592 

347 

>0.02 

0.52813“ 

14.4 ± 1.8 

2.1 ± 1.8 

14.172 

23.8 

106134-3420 

SWASP 

3232-4573 

4700 

>0.03 

12.7 ±2.2 

13.958 

106.4 

J07514-0925 

ASAS 

2623-5131 

198 

>0.1 

0.17832“ 

80.3 ± 20.3 

43.5 ± 20.8 

14.126 

205.9 


Catalina 

3466-6368 

119 

>0.1 

0.17832“ 

12.0 ± 3.2 

3.6 ±3.0 

14.168 

18.7 

J0805-1058 

ASAS 

1868-5168 

570 

>0.1 

0.17370“ 

19.9 ± 3.8 

11.7 ±4.8 

12.270 

65.6 


NSVS 

1488-1630 

132 

>0.01 

0.17370“ 

33.0 ±7.4 

10.9 ±7.2 

12.812 

59.8 

J1356-4934 

ASAS 

1900-5088 

729 

>0.04 

3.3 ±3.9 

12.269 

74.2 

J1411-3053 

ASAS 

1902-5088 

1060 

0.02449 “ 

46.8 ± 3.8 

12.342 

88.2 


SWASP 

3860-4614 

13079 

0.02449 “ 

51.2 ±2.0 

12.723 

165.6 

J16324-0759 

ASAS 

2175-5106 

399 

>0.1 

2.9515“ 

23.9 ±5.1 

3.9 ±5.2 

12.763 

74.5 


Catalina 

3466-6471 

338 

>0.02 

2.9515“ 

49.8 ± 8.4 

7.5 ±9.5 

12.587 

121.6 


NSVS 

1275-1417 

115 

>0.01 

2.9515“ 

37.6 ±8.6 

17.9 ± 8.9 

13.248 

69.3 

J17314-0647 

ASAS 

2727-5009 

73 

>0.1 

1.17334“ 

98.7 ± 49.6 
44.5 ±51.5 

13.799 

299.1 


Catalina 

3466-6457 

105 

>0.1 

1.17334“ 

74.9 ± 4.7 

21.4 ±9.2 

13.825 

67.0 

J17364-2806 

SWASP 

3128-4325 

9140 

1.33320'' 

10.2 ±0.2 

11.639 

17.0 

J1753-5007 

ASAS 

1947-5137 

544 

>0.1 

18.6 ±5.7 

12.955 

94.2 

J1903-3528 

SWASP 

3860-4551 

7388 

>0.01 

3.6 ±0.6 

13.089 

35.2 

J2038-2657 

NSVS 

1348-1483 

42 

1.860'' 

56.4 ± 8.4 

11.950 

51.1 


SWASP 

3958-4614 

8332 

1.87022'' 

12.9 ±0.4 

11.856 

28.2 

J2205-3141 

ASAS 

1873-5166 

521 

0.34156'' 

40.0 ± 5.0 

12.381 

81.4 


Catalina 

3598-6217 

252 

0.34156'' 

46.2 ± 6.2 

12.07 

64.4 


SWASP 

3862-4614 

22731 

0.34156'' 

26.7 ± 1.0 

12.409 

110.1 

J2254-5515 

ASAS 

1869-5168 

672 

>0.02 

1.22702“ 

20.0 ± 3.8 

4.7 ±3.9 

12.113 

70.9 


Catalina 

3580-6076 

78 

>0.1 

1.22702“ 

98.4 ± 16.8 

28.1 ± 15.3 

12.442 

103.1 

J23494-3844 

NSVS 

1321-1579 

261 

>0.1 
0.46252 “ 

19.9 ± 3.9 

5.0 ±3.9 

12.287 

44.7 


SWASP 

3154-4669 

12175 

>0.01 
0.46252 “ 

1.7 ±0.2 

1.1 ±0.2 

11.640 

15.6 


^ Spectroscopic period. 

^ Possible spectroscopic period. 

^ Ellipsoidal variations at half-spectroscopic period. 
Photometric period. 


GALEX J2038-2657 also implies the presence of surface inho¬ 
mogeneities (see Section 3.3.1). 

An analysis of time series helps constrain the nature of the 
companion ( [Maxted et al.|20d^ . For example, a simple geometric 
model suggests that the presence of a late-type companion gener¬ 
ally leads to detectable photometric variations phased on the orbital 
period. This reflection effect scales as R^, where R 2 is the se condary 
radius, but only as where a is the orbital separation jMaxted 
|et al.|2002^ . Interestingly, an application of these simple relations 


confirms the results of detailed light curve modelling and most im¬ 
portantly that, for a given mass function, the effect of a lower bi¬ 
nary inclination, which reduces the light contrast between inferior 
and superior conjunctions as well as its intensity through increased 
binary separation, is compensated by the increased mass and radius 
of the secondary calculated using the mass function, hence increas¬ 
ing the fraction of intercepted light. Following |Maxted et al.|{2002) l 
and adding a slight modification to account for the effect of incli- 
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nation on the visibility of the exposed hemisphere at inferior and 
superior conjunctions, the amplitude of the variations is given by: 


6 m = 2.5 



where 


= 1 + 


Ri \'^^l±sini 

^ V2a' 2 


where /+ is the relative flux at superior conjunction, /“ is at in¬ 
ferior conjunction, i is the binary inclination and Ri is the pri¬ 
mary radius estimated from the measured surface gravity and as¬ 
sumed mass (0.23 or 0.47 Mq). The mass of the putative late-type 
secondary was estimated using the binary mass function, the as¬ 
sumed primary mass and by varying the inclination angle: The ra¬ 
dius is then estimated following the mass-radius relation for late- 
type stars of |Caillault & Patterson| ( [T990| l. Applications of this ap¬ 
proximate formula for 6 m lead to an overestimation of the ampli¬ 
tude of a factor of x: 3 when applied to the well known case of 
GALEX 103214-4727 (Table[^: The semi-amplitude of the phased 
light curve of GALEX J0321-1-4727 reaches 44 and 61 mmag in 
the SWASP and NSVS data sets, respectively, and reveals the pres¬ 
ence of an irradiated late-type companion. This simple model also 
shows that the amplitude of the variations is more or less constant 
(±30 per cent) when varying the inclination as shown in the de¬ 
tailed models of |Maxted et al.|f2002^ . Applying a factor of 0.3 to 
the amplitude calculated with the simple formula for 6 m presented 
above should allow us to confirm or rule out the presence of a late- 
type companion in the new binaries. For example, the companion 
to GALEX J2349±3844 ( [Kawka et al.|2010) is almost certainly a 
white dwarf: The predicted semi-amplitude of variations due to a 
late-type companion is x: 70 mmag, while the observed variations 
are less than 1 and 5 mmag in the SWASP and NSVS phased light 
curves, respectively (Table|^. Based on this insight, the photomet¬ 
ric times series will help constrain in the following Sections the 
nature of the companion in the new binaries. 


3.1.3 Overview of the radial velocity data set 

We measured the radial velocities by fitting a Gaussian function to 
the deep and narrow Her core for most red spectra, or Held 6678.15 
in a few instances described below. In the blue we used the H/3 
core, Heii/14685.698, or Hei/i4471.48 if necessary (see below). All 
measured velocities are heliocentric corrected and tabulated in Ta¬ 
ble |bT] in Appendix B. For each target Table [bT] also includes the 
number of spectra, the average velocity (i;) and velocity dispersion 

(o-O- 

Fig. shows the distribution of the measured velocity disper¬ 
sion. The sample includes three objects published earlier jKawkaj 
jet al.|20l0l|Vennes et al.|2012^ and seven new identifications de¬ 
scribed in the following Section. Three additional objects show 
significant radial velocity variations (cr„ > lOkms^*), but with 
insufficient sampling to determine the orbital parameters. Adding 
two likely close binaries identified through photometric variations 
(GA T. EX J1736+2806 and J2038-2657, see Section 3.1.2) but for 
which we only dispose of radial velocity measurements of the sec¬ 
ondary, we estimate that 15 out of 41 hot subdwarfs presently in¬ 
vestigated are in close binaries, or a 37 per cent yield, lower than 
previously estimated (e.g., jCopperwheat et al.|2011| l. Our survey 
strategy aimed at short-period binaries would be insensitive to long- 
period, low-amplitude variation (< 10 km s“') systems. A detailed 
comparison with the sample of known hot subdwarf binaries should 



Figure 6. Number of objects per velocity dispersion bin (width= 5 kms^*). 
The sample includes data f rom Tab lejBljin Appendix B and published re¬ 
sults from the same survey [Kawka et al.|2010||Vennes et al.|20T^. 


allow us to secure a global estimate of binarity in this population 
(Section 4.1). 


3.2 Individual properties 

Section 3.2.1 describes objects with variable radial velocities sug¬ 
gesting the presence of a close binary companion: Fig. 0Fig.0 
Fig. 1^ Fig. [T^ Fig. fTT) Fig. [T^ and Fig. [T3] show results of the 
period analysis {Ijx^ versus frequency) for this group of objects. 
The confidence level is set at Icr (66 per cent) for a four-parameter 
{p = 4) analysis with the normalized on the best-fitting solution 
and the radial velocity measurements phased on the best-fitting pe¬ 
riod. Sections 3.2.2, 3.2.3, and 3.2.4 review the properties of the 
remaining systems, i.e., those with composite spectra, unresolved 
radial velocity variations, or photometric variability, respectively, 
and the likelihood that they might belong to the close binary popu¬ 
lation. Finally, Section 3.3 presents known facts concerning the re¬ 
maining objects. In the following Section, the subscript “1” refers 
to the hot subdwarf and the subscript “2” refers to its companion. 
Similarly, the suffix “B” designates the companion. The binary pa¬ 
rameters are listed in Table [^including the number of spectra per 
object (A) and the dispersion in velocity residuals (cr,,,.). 


3.2.1 Close binaries 

The sdB star GALEX J0507±0348 is part of a newly identi¬ 
fied spectroscopic binary. The star is close to the ZAEHB and 
may be a low-mass sdB star. The Ha radial velocity measure¬ 
ments are phased on a period of ~0.528d (Fig. |^. The SED of 
GALEX J0507±0348 shows an infrared excess but in the WISE 
W3 band only. The nearby object (sep. = 17 arcsec) visible on 
photographic plates is not likely to affect the WIS E measurements. 
Also, low-dispersion spectra show weak CaH&K lines with an 
equivalent width of £'.VL(CaK)= 270 mA. The calcium doublet 
could indicate the presence of a late-type companion. However, 
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Table 3. Spectroscopic binary parameters. 


Parameter 

J0507+0348 

J0751+0925 

J0805-1058 

J1632+0759 

J1731+0647 

12205-3141 

12254-5515 

p(d) 

0.528127 

0.178319 

0.173703 

2.9515“ 

1.17334 

0.341543 

1.22702 

o-/>(d) 

0.000013 

0.000005 

0.000002 

0.0006 

0.00004 

0.000008 

0.00005 

To (HJD) 

2456315.349 

2455972.827 

2456299.0335 

2456150.701 

2456313.119 

2456313.3387 

2456444.616 

(7 (To) 

0.015 

0.001 

0.0026 

0.016 

0.004 

0.0005 

0.001 

K (kms^*) 

68.2 ±2.5 

147.7 ± 2.2 

29.2 ± 1.3 

54.9 ± 4.6 

87.7 ±4.1 

47.8 ±2.2 

79.7 ± 2.6 

y (kms“*) 

96.2 ± 1.8 

15.5 ± 1.6 

58.2 ±0.9 

-31.6 ±2.7 

-39.1 ±3.0 

-19.4 ± 1.7 

4.2 ± 2.0 

/sec (Mq) 

0.017 ± 0.002 

0.059 ± 0.003 

(4.4 ± 0.6) X lO^'* 

0.048 ±0.013 

0.080 ± 0.012 

0.0037 ± 0.0005 

0.063 ± 0.006 

Ml (Mq) 

(0.47) 

(0.47) 

(0.234) 

(0.47) 

(0.47) 

(0.47) 

(0.47) 

M 2 (Mq) 

> 0.20 

> 0.34 

> 0.03 

>0.31 

> 0.39 

> 0.11 

> 0.35 

N 

16 

19 

23 

12 

16 

13 

24 

(Tp (km s“^) 

6.5 

6.4 

2.8 

8.6 

6.8 

5.4 

8.2 

Notes 

probable WD 

probable WD 

low mass sdB, 

probable 

probable 

reflection, 

probable 


secondary 

secondary 

possible BD 
secondary 

WD secondary 

WD secondary 

dM secondary 

WD secondary 


Based in part on the period obtained by 


Barlow et al. 


2014 


:P = 2.951 ±0.001 d. 


GALEX J 0507±0348 


GALEX J0751±0925 





Frequency (d^*) 


Frequency (d^') 


Figure 7. (Bottom) period analysis of radial velocity measurements of 
GALEX J0507±0348 (full line) and 66 per cent confidence level (dashed 
line). (Top) radial velocity measurements phased on the orbital period 
(0.52813 d) and and best-fitting sine curve (full line) with the dispersion in 
velocity residuals shown in upper-right. The initial epoch Tq corresponds to 
inferior conjunction of the sdB star. Details are presented in Section 3.2.1. 


no radial velocity measurements were obtained in that spectral re¬ 
gion and we could not confirm variations in the line position. More¬ 
over, we could not confirm the presence of other late-type spectral 
signatures such as Mgi lines, and our composite spectral analysis 
( [Nemeth, Kawka & Vennes|2012| > rejects the presence of a com¬ 
panion with a flux contribution above 1 percent in the optical range. 
A series of high dispersion spectra is required to determine whether 
the CaH&K lines originate from the companion, the interstellar 
medium (ISM), or in the circumstellar environment. The mass func¬ 
tion allows us to infer M 2 > 0.20 Mq assuming Mj = 0.47 Mq, or 
M 2 > 0.15 Mo assuming M| = 0.30 Mq. A late-type (dM,dK) com¬ 


Figure 8 . Same as Fig.^but for GALEX J0751±0925. 

panion would satisfy these constraints. However, the Catalina time 
series folded on the orbital period constrains the photometric vari¬ 
ations to a semi-amplitude lower than 2 mmag. Reflection effect on 
a late-type star with a mass exceeding 0.2 Mq would result in vari¬ 
ations of a semi-amplitude of ~ 60 mmag as observed in the case 
of GALEX J0321+4727 jKawka et al.|M^ . We conclude that the 
companion is most likely a white dwarf: At a binary inclination 
i < 30° the mass function implies a minimum mass of 0.51 Mq that 
would be consistent with a normal white dwarf star. 

The sdB star GALEX J0751+0925 is part of a close binary 
with the largest velocity semi-amplitude measured in the present 
sample, K ~148kms“* (Fig. [^. The mass function implies the 
presence of a companion relatively more massive than in other sim¬ 
ilar systems with M 2 > 0.34 Mq assuming M| = 0.47 Mq. The 
SED of this system also appears to show an infrared excess in the 
Wl, W2, and W3 bands which may be caused, in part, by a nearby 
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Figure 9. Same as Fig.^but for GALEX J0805-1058. 


GALEX J1632 + 0759 




Frequency (d"^) 


GALEX J1731+0647 



Figure 11. Same as Fig. |7]but for GALEX J 1731+0647. 


GALEX J2205-3141 



Figure 10. Same as Fig. [Tjbut for GALEX J1632+0759. The tick mark 
above the best-htting period indicates the results of the period analysis of 
[Barlow et ahlpoTi) , P = 2.951 ± 0.001 d. 

star only 6 arcsec away (particularly in the W3 band). The Catalina 
and ASAS light curves do not show significant variations at the 
orbital period. Again, a comparison with the photometrically vari¬ 
able system GALEX J0321-)-4727 indicates that the companion is 
probably a white dwarf. The semi-amplitude of the phased light 
curve of GALEX J0751+0925 is limited to 4 mmag in the Catalina 
observations although the minimum secondary mass in this shorter- 
period system is larger than in GALEX J0321+4727, i.e., 0.34 ver¬ 
sus 0.13 Mq, and the predicted semi-amplitude of variations due to 
a late-type companion would be x: 190 mmag. Therefore, the ab¬ 


Figure 12. Same as Fig. [Tjbut for GALEX J2205-3141. The photometric 
period is marked close to the peak frequency of the velocity periodogram. 

sence of a reflection effect in GALEX J0751+0925 rules out the 
presence of a late-type companion leaving only the possibility of a 
> 0.34 Mq white dwarf companion, or > 0.50 Mq if i < 50°. 

The sdB GALEX J0805-1058 clearly lies below the ZAEHB 
(Fig. 1^; following the evolutionary tracks of |Driebe et al.|p998| l 
the mass of the subdwarf is estimated to be in the range 0.2-0.3 Mq. 
The low velocity amplitude (Fig. [^ and small mass function im¬ 
ply a very low mass for the companion, M 2 > 0.03 Mq, assum¬ 
ing Ml = 0.234 Mq. At inclinations higher than i = 26°, the sec¬ 
ondary mass remains lower than 0.08 Mq and the object is sub- 
stellar. Assuming a probability distribution for the inclination an- 
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GALEX J2254-5515 



Figure 13. Same as Fig. |7]but for GALEX J2254-5515. 

gle i of the form P, di = sin i di, inclinations higher than 26° have 
P{> 26) =89.9 per cent probability of occurring. At lower inclina¬ 
tions (10° < i < 26°, i.e., P(10 - 26) =8.6 per cent), the secondary 
mass does not exceed 0.3 Mq and the object would be a low-mass M 
dwarf. At very low inclinations (i < 7°, i.e., P(< 7) =0.7 per cent) 
the secondary would be a normal white dwarf star (M > 0.5 Mq). 
The SED shows a single, hot subdwarf star, i.e., as far as the WIS E 
W2 band, and the faint (Am ~6 mag), nearby (sep. ~11 arcsec) ob¬ 
ject visible in photographic plates does not appear to contaminate 
the SED. The ASAS and NSVS light curves do not show signifi¬ 
cant variations when folded on the orbital period, i.e., < 12 mmag, 
while the predicted semi-amplitude of variations due to a substellar 
object {R 2 ~ 0.1 Pq, see a review by |Chabrier et al.|2009| l would 
be 20 mmag. Variations of the order of 10 mmag have been ob¬ 
served in brown dwarf plus hot subdwarf binaries (see, e.g.,|S ch^ 
[fenroth et al.||2014^ and such variations would be detectable in 
GALEX J0805-1058 in quality photometric time series. We con¬ 
clude that the apparent lack of variations is a consequence of the 
small radius of a substellar secondary. 

[Barlow et al.| ( [2012^ recorded radial velocity variations in 
spectra of the sdB star GALEX J1632-I-0759. Their data suggested 
a period ranging from 2 to 11 days. Our measurements, based on 
Ho- and Hei5875.621 in the red, and H/3 and Hei4685.698 in the 
blue, also revealed large velocity variations (Fig. [10^ . Recently, 
[Barlow et al.j j2014^ obtained new radial velocity measurements 
and determined a period of 2.951 ± 0.001 d. We restricted our pe¬ 
riod analysis to frequencies between 0.3 and 0.4 d“* and recovered 
an identical period. The mass function implies a secondary mass 
M 2 > 0.31 Mq, assuming Mi = 0.47 Mq. The SED shows a large 
flux excess apparent in the 2MASS and WISE bands as well as 
heavy extinction in the ultraviolet range. The measured extinction 
coefficient (Eg-v = 0.4) largely exceeds the coefficient inferred 
from the maps of |Schlegel, Finkbeiner, & Davis[ ( [1998^ , Eg^y = 
0.08. The additional extinction probably originates in the imme¬ 
diate, possibly dusty, circumstellar environment of the system. An 
inspection of our acquisition images of GALEX J1632-I-0759 re¬ 
veals the presence of a nearby star; we measured a separation of 


2.3 arcsec at a position angle of 225°. The 2MASS and WISE 
photometric measurements are likely contaminated by this object. 
[0stensen et al.[ ( [2005[ l also resolved GALEX J1632-I-0759 and the 
nearby star and measured a separation of 2.1 arcsec. In addition 
to GALEX J1632-f0759, [Barlow et al.|j2014^ obtained radial ve¬ 
locity measurements of the nearby star which they classified as 
a late G dwarf or early K dwarf that is itself in a close binary: 
The radial velocity varied with a period of 1.42 ± 0.01 d. |Barlow| 
[et al.[ j2014[ > also found that the systems share the same systemic 
velocity suggesting that this is a quadruple system. The ASAS, 
Catalina, and NSVS time series constrain photometric variations 
to semi-amplitudes lower than 4, 8, and 18 mmag. The predicted 
semi-amplitude of photometric variations due to the presence of a 
late-type companion would be 60 mmag. We conclude that the 
secondary star is most probably a white dwarf with a mass ranging 
from 1.3 to 0.5 Mq assuming a low inclination (24 < i < 46°), or 
with a peculiar low mass (0.3 - 0.5 Mq) assuming i > 46°. 

The new binary system GALEX J1731-I-0647 (Fig. pT) har¬ 
bours the heaviest binary companion identified in our sample. The 
mass function implies a mass M 2 > 0.39 Mq, assuming Mi = 
0.47 Mq. The field surrounding this suhdwarf is relatively crowded 
but only two objects are found between 13 and 15 arcsec away and 
with photographic magnitude differentials Am ~3 and 5 mag. These 
objects would not affect the SED which shows a single hot subd¬ 
warf. The lack of photometric variations, < 45 mmag in ASAS 
time series and <21 mmag in Catalina time series, compared to 
expected variations of « 90 mmag due to a relatively large M or 
K dwarf suggests that the companion is most likely a white dwarf. 
We infer a mass between 1.3 and 0.5 Mq assuming an inclination 
of 29 < i < 58°, or a peculiar low mass (0.4 - 0.5 Mq) assuming 
i> 58°. 

GALEX 12205-3141 is a close binary with P x 0.34 d 
(Fig. [12^ showing a reflection effect in the SWASP time series 
(semi-amplitude Am a; 27 mmag). Similar variations were also 
observed in the ASAS and Catalina time series. The mass func¬ 
tion is consistent with the presence of a late-type companion 
(M2 > 0.11 Mq). Photometric time series from the Catalina survey, 
SWASP (ISWASP 1220551.98-314103.9), and ASAS show vari¬ 
ability with mutually consistent periods of 0.341559 ± 0.000003, 
0.341563±0.000002, and 0.341561 ±0.000002d, respectively. The 
photometric periods are somewhat longer than the spectroscopic 
orbital period P = 0.341543 ± 0.000008 d. The radial velocity 
measurements are based on the Hei 4471.48 and 6678.15 lines. 
We noted that the Balmer Ha and Hfl lines cores are asymmetric 
and are possibly contaminated by emission from the companion as 
noted in the case of AA Dor ( [Vuckovic et al.|2008[ ). Fig. [14| shows 
the SWASP measurements phased on the photometric period. We 
identify the initial epoch with the passage of the secondary star 
at superior conjunction corresponding to maximum reflected light. 
Although the photometric variations are clearly caused by the re¬ 
flection of the primary light on a late-type dwarf companion, the 
phasing error between photometric and spectroscopic ephemeris is 
AO 0.1. We attribute this error to a large gap between the epoch 
of the spectroscopic observations and that of the photometric ob¬ 
servations. The SED shows a mild flux excess in the IR to mid- 
IR range possibly due to the late-type companion. A star found 
9 arcsec away and 4 mag fainter does not affect the SED. How¬ 
ever, renormalizing on the J band rather than the V band nearly 
eradicates this excess. Assuming a possible K band contribution 
from the companion of 15 to 40 per cent, we estimated for the M 
dwarf companion Mg^^ ~ 7.5 to 6.5 if Mg i x 5.5. Bearing in mind 
that reprocessing of ultraviolet radiation from the hot primary into 
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Figure 14. (Left panels) Fourier transform analysis of the variable stars GALEX J2205-3141 (top), J1736+2806 (middle) and J2038-2657 (bottom), 
and phased light curves (right panels). The identification of the photometric period with the spectroscopic period clearly indicates that the light curve of 
GALEX J2205-3141 shows the effect of reflection of the bright primary on the secondary star. The initial epoch Tq in GALEX J2205-314I comesponds 
to the passage of the secondary star at superior conjunction. Without knowing their orbital periods, we can only state that the cool, secondary stars in both 
GALEX J1736+2806 and J2038-2657 are variable. The photometric periods are marked with star symbols: P = 0.341563 d (12205-3141), 1.333204 d 
(11736+2806), and 1.870221 d (12038-2657). 


the cool secondary atmosphere should contribute to this IR excess, 
the absolute K magnitude of the secondary star corresponds to a 
spectral type later than M3-4, or a mass M 2 < 0.24-0.4 Mq which 
requires an orbital inclination i > 20-30°. We find possible evi¬ 
dence of extinction in the ultraviolet range in excess of the extinc¬ 
tion expected from the |Schlegel, Finkbeiner, & Davis|(T998^ map, 
although the GALEX NUV dip may be the result of larger uncer¬ 


tainties than estimated. This system is the only confirmed binary in 
our sample comprised of a hot subdwarf and late-type companion. 

The sdB GALEX J2254-5515 shows large radial velocity 
variations (Fig. |13| l although the Catalina and ASAS time series 
indicate that the star is not photometrically variable with semi¬ 
amplitudes lower than 28 and 5 mmag, respectively. The minimum 
mass of the secondary, M 2 > 0.35 Mq assuming Mi = 0.47 Mq, 
combined with the lack of photometric variability when compared 
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to expected variations of 150 mmag caused by a reflection effect on 
a putative late-type companion imply that the companion is a white 
dwarf. 

We neglected the possible effect of orbital eccentricity in the 
period analysis. The orbits of post-common envelope binaries is 
expected to be circular due to the synchronization during the post- 
CE phase. However, eccentric orbits in close binaries containing 
a subdwarf were reported by [Edelmann et al.| (|2005f and |Kawka| 
|et al.| ( [2012a[ >. In these cases the eccentricity was small (e < 0.1). 
Larger eccentricities were reported for long period binaries, such 
as BD+20°3070, BD-h 34°1543, Eeige 87 ^Vos et al.||2013f and 
PG 1449+653 jBarlow et al.|2013^ . Eccentric orbits may indicate 
the presence of a circumbinary disc ( [Artymowicz et al.|199H l. 

Now, we summarize additional constraints on the properties 
of spectroscopic composites, other likely systems showing radial 
velocity variations, and systems displaying photometric variability. 

3.2.2 Composite spectra 

Using spectral decomposition, [Nemeth, Kawka & Vennes| (2012^ 
classified GALEX J0047+0337 as a binary consisting of a hot sdB 
and a main-sequence F star. The radial velocity measurements ob¬ 
tained for GALEX J0047+0337B imply a constant velocity with 
standard deviation of only 6.3 km s“* and include a single measure¬ 
ment deviating from the average velocity hy more than lOkms^L 
The AS AS and NS VS photometry do not show evidence of signif¬ 
icant variations: The AS AS data constrain potential variations to 
a semi-amplitude lower than 4.4 mmag for all periods larger than 
0.5 hr. The EFOSC acquisition images revealed a nearhy star ap¬ 
proximately ~3 arcsec away at a position angle of 344°. The object 
is about 1.2 mag fainter in R and the quoted WISE and 2MASS 
magnitudes include both stars since they would not be resolved in 
either surveys. Fortunately, our optical spectra were not contami¬ 
nated by the nearby star and the composite nature of the object is 
not affected. 

GALEX J1411+7037 and J1753-5007 are sdB stars with F- 
type companions. Their SEDs are consistent with the presence of 
a luminous companion derived from the spectral decomposition of 
[Nemeth, Kawka & Vennes| ( |2012[ . The Ha line profile in each star 
is dominated by the main-sequence star and no significant radial 
velocity variations have been found for these objects. The ASAS 
times series of GALEX J1753-5007 constrain photometric varia¬ 
tions to a semi-amplitude lower than 19 mmag. 

The SED of GALEX J1356-4934 shows significant infrared 
excess. An inspection of the acquisition images did not reveal a 
resolvable, nearby companion and radial velocity measurements 
show only marginal variability with radial velocity maxima reach¬ 
ing a span of 20 km s“*. First, we performed a SED decomposition 
to estimate the spectral type of the companion. We adopted the sdB 
parameters determined by [Nemeth, Kawka & VennesH2012) and 
calculated sdB absolute magnitudes of Mk - 5.47 and My = 4.46. 
Adopting the apparent visual magnitude V = 12.3 and 2MASS 
magnitude K = 11.633 and using the main-sequence colour and 
absolute magnitude relations from jPecaut & Mamajek| ( |2013[ >, we 
determined the absolute visual and infrared magnitudes of the late- 
type companion, Mk 2 = 3.57 and Mv .2 = 5.36, and a distance of 
444 pc. Consequently, the companion mass is 0.94 Mq correspond¬ 
ing to a G8V star. Next, we performed a spectral decomposition 
with XTGRID (Nemeth, Kawka & Vennes||2012| > making use of 
both the blue and red spectra of GALEX J1356-4934. The spec¬ 
tral decomposition showed that the companion contributes 27 per 
cent of the flux at 7000 A. The new parameters of the sdB star are 


Teff = 32370!|o K, logg = 5.72!°;?’, logHe/H = -2.75!“;“ and 
do not differ significantly from our earlier measurements. The pa¬ 
rameters of the companion are T^ff = 5470, logg = 4.47, [Fe/H] 
= 0.003, also corresponding to a G8 main-sequence star. These 
values supersede those of [Nemeth, Kawka & Vennes|p01^ for 
GALEX J1356-4934. The ASAS time series limits the photomet¬ 
ric variations to a semi-amplitude of 3 mmag. 

Optical spectra of subdwarf plus early-type F-stars are domi¬ 
nated in the red by the companion. Because the mass ratio is >3, 
high-dispersion spectroscopy is required to detect the secondary 
star motion. 


3.2.3 Radial velocity variable 

Other objects, in addition to the confirmed binaries listed in Tablej^ 
are likely close systems. The measured radial velocity extrema sug¬ 
gest that these subdwarfs are in close orbit with a companion, but 
the small number of spectra did not allow us to perform a period 
analysis. 

We measured velocity extrema At; 80kms“* for 

GALEX J0613+3420, Av ^ 28kms-' for GALEX 10812+1601, 
and Av « 35 km s ' for GALEX 11903-3528. The SWASP time 
series for GALEX 10613+3420 and GALEX 11903-3528 con¬ 
strain photometric variations to maximum semi-amplitudes of 13 
and 4 mmag, respectively, which exclude the presence of close late- 
type companions. Further investigations are required to clarify their 
binary status. 

The SED of each object does not reveal the presence of a com¬ 
panion, however the SED of GALEX 10613+3420 shows evidence 
of a large interstellar extinction (Schlegel, Finkbeiner, & Davis| 
|1998| >, and a possible excess (Eg-v = 0.64) above the interstellar 
value (Eg-v = 0.36). 


3.2.4 Photometrically variable 

The SWASP light curve of the sdB plus F7V pair 
GALEX 11736+2806 (ISWASP 1173651.18+280634.6) varies 
with a period P = 1.33 d and a semi-amplitude of 11 mmag 
(Fig. [14}. No significant variations were observed in a nearby 
comparison object (ISWASP 1173635.80+280902.2). The group¬ 
ing of data points observed in the light curve are also observed 
in the light curve of the nearby object and, therefore, it must be 
an artefact of data sampling. Using the SED we found that the 
absolute V magnitude of the companion is ~0.72 mag brighter 
than the sdB star consistent with a value of ~0.81 mag obtained 
by [Nemeth, Kawka & Vennes| ( |2012} . The absolute magnitude of 
a late F7 star. My ~4, would imply for the hot subdwarf My ~4.7; 
The atmospheric parameters of the hot subdwarf are very uncertain 
(Nemeth, Kawka & Vennes[[2012[ > but would be reconciled with 
the companion spectral type at the lowest acceptable temperature 
(30,000K at logg = 5.7). The photometric variability may be 
caused by irradiation of the exposed hemisphere of the F star 
although we failed to detect radial velocity variations at the same 
period. 

GALEX 12038-2657 is a relatively luminous hot sdO star 
with a G type companion (Nemeth, Kawka & Venn^[2012[ >. 
Our spectroscopic observations revealed variability in the Ha pro¬ 
file (Fig. [T5) on a time-scale of a day or less. However, cross¬ 
correlation measurements in the spectral series dominated by the G 
companion show little variations with a dispersion cr,, = 7.8 km s“' 
comparable to the expected accuracy of the wavelength scale. The 
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measurements imply that the velocity semi-amplitude of the G8III 
star does not exceed x: 16kms“'. Fig. 16 shows the SED of 
GALEX J2038-2657 where the ultraviolet range is dominated by 
the hot subdwarf and the optical range by the red giant. The SWASP 
time series (ISWASP J203850.49-265754.2) reveals variations of 
12 mmag semi-amplitude over a period of 1.87022 d (Pig.|14^ that 
are confirmed by similar variations in a short NSVS time series. 
Again, no significant variations were observed in a nearby com¬ 
parison object (ISWASP J203851.22-265943.1). These variations 
are most likely linked to the observed spectroscopic variability, but 
they cannot yet be clearly associated to a possible orbital period. 

The system shares some properties with the sdB plus KIII- 
IV system HD 185510 ( [Jeffery et al.[T9^|Pekel et al.|1993| > and 
the sdO plus KO III system EE Aqr (Vaccaro & Wilson|2003p. All 


three systems have an evolved secondary star, from sub-giant to 
giant, and all three are photometrically variable. However, the hot 
subdwarf in HD 185510 is possibly the progenitor of a low-mass, 
helium white dwarf (0.3 Mq, [Jeffery & Simon|1997^ . Photomet¬ 
ric variations in HD 185510 and EE Aqr coincide with the orbital 
period and are caused by irradiation of the exposed hemisphere of 
the secondary stars. Moreover, the orbital periods HD 185510 and 
EF Aqr are 20.7 and 9.2 d, respectively, with orbital separations 
of ~43 and ~25 Rq, respectively, and well outside the radius of a 
sub-giant or giant star (R(K0 III) 16 Rq). 

Without an estimate of the orbital period we can only set lim¬ 
its to the orbital parameters, such as the binary separation. The 
identification of the late-type giant secondary is based on spectral 
decomposition ( [Nemeth, Kawka & Vennesj2012^ : The absolute V 
magnitude of the hot sdO star is only about ~ 1.0 mag fainter than 
its companion. Adopting a G8 III type from the spectral decompo¬ 
sition shown in Fig.[^ the absolute magnitude of the companion is 
Mv/(G8III)= 0.9, implying an absolute magnitude Mv(sdO)= 1.9 
for the primary in agreement with the estimate of [Nemeth, Kawk^ 
[& Venn^ p012[ l, Afv/(sdO)= 2.0. The minimum orbital period 
for a systemic mass of 2-3 Mq and an orbit outside the G8 III ra¬ 
dius (15 Rq) is P >3.1 d. Adopting a radius of 15 Rofor the G8 III 
star, the photometric period of 1.87 d implies a rotation velocity 
Urot ~ 350kms“*. The narrowest features in the SSO spectra have 
a width of t),.(,tsini = 130kms“', that would enforce a low incli¬ 
nation i < 22°. High dispersion spectroscopy is necessary to help 
determine the orbital parameters and help clarify the origin of the 
photometric variations. The most likely scenario is that the photo¬ 
metric variations are caused by a surface spot coupled to the rota¬ 
tion of the star, and that the orbital period probably exceeds several 
days with a low velocity amplitude {K <20kms“*). 


3.3 Notes on other objects from this survey 

GALEX J0047+0958 (HD 4539) is a well known hot sdB star (see, 
e.g., |Kilkenny|1984|l. Spectropolarimetric measurements hint at the 
presence of a weak magnetic field (~ 0.5 kG [Landstreet et al.[2012[ l. 
[Schoenaers & Lynas-Gray] ( [2007[ reported line profile variations 
and radial velocity variations of a few kms“* that may be due to 
g-mode pulsations. [Lynas-Gray[ ( |20T^ obtained photometric series 
and measured variations with a frequency of 9.285 ± 0.003 d“' and 
an amplitude of 0.0023 ± 0.0003 mag. This photometric frequency 
is consistent with one of the frequencies (9.2875 ± 0.0003 d“*) de¬ 
termined from low-amplitude radial velocity variations, and both 
are possibly associated to stellar pulsation. 

The SED of GALEX J0049-H2056 (Eig.|^ shows an IR ex¬ 
cess that could be attributed to a yet unidentified companion or 
nearby object, or to a dusty environment. 



x(A) 


Figure 15. Spectra of the photometrically variable sdO+G8III star 
GALEX J2038-2657 obtained at SSO and La Silla and showing short-term 
variable Hu emission. 



a(A) 


Figure 16. Spectral energy distribution of GALEX J2038-2657 combin¬ 
ing a cool G8III secondary and a sdO primary star (full line). Individual 
contributions are shown with grey lines. The effect of interstellar extinction 
(Eb-v = 0.08) is included. 


The sdB GALEX J0059+1544 (PHL 932) is embedded in 
an emission nebula. However, [Erew et al.[[20ldl l have shown that 
the association is only coincidental, but that PHL 932 does con¬ 
tribute and ionize a dense region of the ISM surrounding it. The 
SED of this object shows, as in the case of GALEX J0049-(-2056, 
a considerable IR excess (Fig. [Al[ (. Several radial velocity mea¬ 
surements of PHL 932 were reported in literature. [Arp & Scargle[ 
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1 1967^ measured 15 ± 20kms * using two low-dispersion spec¬ 


tra. 


Edelmann 


j2003^ measured 18±2kms“' using echelle spectra. 
These velocities are in agreement with our measurements (D = 16.7, 
cr = 3.1 km s“*) and, therefore, it does not appear that PHL 932 is in 
a close binary.|Geier & Heber|(|2012[) report a rotational velocity of 


Landstreet et al. 


2012|l obtained spectropo- 


Vrot = 9.0 ± 1.3 km s“*. 
larimetric measurements of PHL 932 but did not detect a magnetic 
field with an upper limit of ~300 G. 

[Brown et 'ei] ( |2008l classified GA l. F.X J0206-H1438 
(CHSS 3497) as a hot subdwarf. Our radial velocity measurements 
vary only marginally (cr„ < lOkms^'), and we do not dispose of 
sufficient data to determine a period. A radial velocity measure¬ 
ment of Vr = 7 ± 16 km s“* was obtained by 1 


Brown et al. 


(2008; 


and 


is consistent with our measurements (v = 13.8, cr„ = 7.5 km s '). 


GAl.F.X J0232-H4411 (FBS 0229-H439) was classified as a sdB 
star in the First Byurakan Survey of blue stellar objects jMickaelian| 
|2008| l. 

jCopperwheat et ar] ( |2011^ presented a set of radial velocity 
measurements for GALFX J0401-3223 which suggest that the sdB 
star is in a close binary system, but were unable to determine the or¬ 
bital period with limited data. |Copperwheat et al.| ( [2Qll^ measured 
an average velocity and dispersion of u ± cr„ = 55.2 ± 4.4kms“*, 
consistent with our own measurements. We have combined the 
jCopperwheat et aL] ( |2011) data with ours and conducted a period 
search. We found a best period of 1.8574 d, however two signifi¬ 
cant aliases at P = 0.64 and 0.066 d cannot be ruled out. The veloc¬ 
ity semi-amplitude at all three periods does not exceed 10 km s“' 
and excludes a white dwarf or late-type companion. The relatively 
short period and low velocity amplitude imply a minimum mass in 
the substellar range, 0.01-0.04 Mo. The SWASP data folded on the 
best period (1.8574 d) constrain photometric variations to a semi¬ 
amplitude of only 1 mmag, or 8 mmag when folded on any periods 
larger than 0.01 d. The expected variations due to a substellar com¬ 
panion would be as low as 6 mmag at the two longest periods or 
20 mmag at the shortest, but are all significantly larger than the 
SWASP limit. It is not possible to describe the companion with 
present data, although a substellar companion is a distinct possibil¬ 
ity- _ _ 

|0stensen et al.j j2010a| l obtained series of photometric obser¬ 
vations of GALFX J0500-I-0912 in order to search for pulsations 
and concluded that it is not photometrically variable. Our limited 
radial velocity data set does not indicate variability. 

An inspection of the acquisition images of 
GALFX J0657-7324 shows a nearby companion and, therefore, 
the 2MASS and WISE colours of the hot subdwarf are certainly 
contaminated. jHeintzj reported that GALFX J0657-7324 

(HFI 714) is a visual double star with a separation of 1.9 arcsec, 
and our own acquisition image locates the companion 1.8 arcsec 
away at a position angle of 270°. Also, our optical spectra do 
not appear to be contaminated by this object and do not indicate 
variability. 

GALFX J1845-4138 is a relatively cool He-rich subdwarf 
displaying a strong Hei line series and weaker Balmer lines. 
The velocity measurements based on Hei6678.154 (v = -59.7 ± 
3.3 kms^*) are consistent with the measurements based on Ha and 
do not suggest any variability. 

GALFX J1902-5130 is a helium sdO star. [Landstreet et al.j 
( [2012^ obtained spectropolarimetry of GALFX J1902-5130 with a 
measurement that shows that this star does not have a magnetic field 
down to a few hundred gauss. Our radial velocity measurements 
suggest there may be long period, low-amplitude variations. The 


measurements are based on Hen 6560.088A. The object is very 
hot and our spectra display He i emission. 

GALFX J1911 -1406 is also a very hot He-rich subdwarf. The 
velocity measurements are based on He ii6560.088A. 

[Geier & Heberj ([2012J report a rotational velocity of Vrot = 
8.6±L8kms~‘ for GALFX J2153-7003. |Copperwheatet al.| l |2011 1 
obtained several radial velocity measurements of this star and found 
that it is not variable. Their average velocity and dispersion, 39.4 ± 
7.5kms“', are in a close agreement with our own measurements 

(43.4±4.2kms-‘). _ 

GALFX J2344-3426 is a well known sdB star. jOstensen et al.| 
([2010a[> obtained photometric series of the star and found it to be 


Mathys et al. 

2012 

1 and 

Landstreet et al. 


of Lot = 7.3 ± l.Okms^'. 

\2Q12) obtained spectropolarimetry of the star and constrained the 
longitudinal field to 261 G and 246 ± 232 G, respectively. 


4 DISCUSSION 

The new binary identifications are placed into context with a com¬ 
pilation of all known spectroscopically identified binaries (Ap¬ 
pendix C). Table [CT] lists the orbital parameters of these systems. 
The compilation includes hot subdwarfs with an unseen companion 
and spectroscopically identified late- to early-type companions in a 
range of periods from 0.05 to 1363 d. Table [C2| lists the properties 
of the primary as well as kinematical properties of the systems. 

Throughout this discussion we assume for most subdwarfs a 
mass of 0.47 Mq with a few exceptions such as the ELM progeni¬ 
tors (0.23 Mq) and hot sdO companions to early-type stars (1 Mq). 
Using pulsating properties of sdB stars and binary systems for 
which the sdB mass was measured [Fontaine et al.j {2Q\2) found 
that the average mass of a sdB star is 0.47 Mq with a standard de¬ 
viation of 0.031 Mq. [Zhang et al.[ ( [2009[ l found that most sdB stars 
have a mass between 0.42 and 0.54 Mq and an average mass of 
about 0.50Mq. [Zhang et al.jpOO^ used evolutionary models and 
the parameters of a sample of 164 sdB stars. 


4.1 Properties of known binaries: Period and mass function 

Fig.[T^(top) shows the cumulative distribution of orbital periods in 
the population of binaries with a hot subdwarf primary. The deriva¬ 
tive of the function with respect to the logarithm of the period pro¬ 
vides an estimate of the period distribution (Fig. \n\ bottom). Sev¬ 
eral peaks stand-out, particularly at 0.1, 0.5-2.0, 10, and 1000 days. 
The last two peaks are clearly separated by a gap within which few 
binaries are known: A few Be stars with a hot subdwarf compan¬ 
ion populate the gap and the distribution includes spectroscopically 
(UV) confirmed Be-l-sdO ( [Peters et al.[[2008| [2013[ l. We noticed 
a hint of a hierarchy in the period distribution: The shortest peri¬ 
ods coincide with dM companions emerging from a CE phase (low 
mass ratio M 2 /M 1 < 1/2), white dwarfs (M 2 /Af| = 1), and, at high 
mass ratio (M 2 /M 1 > 2), subdwarfs with a subgiant/giant or Be 
companion, and, finally, subdwarfs with a FGK companion at the 
longest periods and emerging from a RLOF. Following a RLOF, 
the orbital separation increases the least for more massive compan¬ 
ions. It is remarkable that main-sequence A-type star companions 
are still missing although they are predicted in population syntheses 
( [Han et al.|2003[ (. Subdwarfs with A-type main sequence compan¬ 
ions should be detectable as UV excess objects or as low-amplitude 
radial velocity variations similar to Be-(-subdwarf binaries. In sum¬ 
mary, binaries near the main peak are mostly white dwarfs plus 
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Figure 17. Cumulative function of period (top), A^< versus log P, and its first 
derivative (bottom). The derivative was smoothed with a Gaussian function 
{FWHM = 0.1 dex), and shown in the upper right comer. 



p (d) 

Figure 18. Measured velocity amplitude versus period with a sub-sample 
of eclipsing binaries shown with open circles. Full lines are labelled with 
the mass of secondary stars which were computed for i = 90°. 

subdwarfs after possibly two episodes of mass-transfer. Note that 
the orbital parameters of many systems with large velocity ampli¬ 
tudes remain unresolved (see, e.g., |Copperw heat et al.|2011[|Geier| 
|et al.|201 la^ and are not included in this analysis. 

Fig.|18|shows the sample of known binaries in the velocity am¬ 
plitude versus period plane. Most eclipsing systems have secondary 
masses between 0.08 to 0.15 Mg. Systems with known white dwarf 
secondaries have secondary masses close to or above 0.60 Mg. 



P (d) 

Figure 19. Same as Fig. |18| but with reflection binaries shown with open cir¬ 
cles. The velocity scale is corrected for the effect of an average inclination 
of 57°. Secondary masses cluster between 0.08 and 0.30 Mg. 

Fig. [T^ also shows the sample of known binaries in the veloc¬ 
ity amplitude versus period plane but with the velocity scale cor¬ 
rected for an average inclination of 57°. The correction allows to 
draw class properties but should not be applied to individual ob¬ 
jects. Secondary masses for systems showing a reflection effect 
range, with the exception of FF Aqr and HD 185510, from 0.08 
to 0.30 Mg. Remarkably, secondary masses for most non-reflecting 
systems cluster near 0.60 Mg and the unseen objects are probably 
white dwarfs. Secondary stars in the long-period range and with 
masses in excess of 0.60 Mg are identifiable as G and K stars. All 
eclipsing systems with K < lOOkms^*, i.e., with an estimated 
M 2 < 0.3 Mg, also show a reflection effect indicative of a late- 
type secondary, while the remaining systems cluster at a higher 
secondary mass M 2 ~ 0.6 Mg and almost certainly harbour a white 
dwarf secondary. Fig.|20|shows secondary mass distribution assum¬ 
ing average system inclination of 57°. This distribution may be de¬ 
scribed by a superposition of two power laws: A shallow distribu¬ 
tion with M 2 > 0.08 Mg, i.e., a = 1.3 between 0.08 and 0.5 Mg 
and a = 2.3 above 0.5 Mg following the initial mass function 
((m) oc rtT'^ of |Kroup^ ( |2001^ , and a steeper distribution (a x 6) 
with M 2 > 0.48 Mg. The former power law encompasses mostly 
M-type dwarfs, many of them showing a reflection effect, while 
the latter encompasses white dwarfs in the 0.5-1.0 Mg mass range. 
This simplified white dwarf distribution represents well the peak 
and high-mass tail of the white dwarf mass distribution (see, e.g., 
[Kepler et al.||2007) but excludes possible low-mass white dwarfs 
(< 0.48 Mg). On the other hand, the late-type stellar mass distri¬ 
bution follows the initial mass function and the expectation of a 
randomly drawn set of late-type stars. The reflection effect is com¬ 
mon in short-period binaries (P < 0.5 d) but is relatively rare at 
longer periods (see Fig.|19|l due to increased binary separations and 
weaker photometric variations: The actual late-type mass distribu¬ 
tion appears as a scaled-up version of the secondary mass distri¬ 
bution in the sub-sample of reflection binaries, but it also includes 
longer period binaries with an indiscernible reflection effect. Note 
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Figure 20. Mass distribution of all known binaries with a hot subdwaif 
primary star as a function of the secondary mass, assuming an average in¬ 
clination of 57° (full histogram). The peak distribution of low-mass stars 
is marked “dM” and that of white dwarfs, “WD”. Binaries showing reflec¬ 
tion elfect in their light curves are shown with a dashed histogram. The full 
lines show synthetic distributions smoothed to two-bins width for a com¬ 
bination of late-type stars and white dwarfs (double-peaked full line), and 
that excluding white dwarfs (dashed line). 


that a third narrow peak is possibly present at « 0.3 Mq. Since this 
peak is mostly made up of companions that do not show the reflec¬ 
tion elfect, the origin of this peak maybe be low-mass white dwarfs. 
It may also be due to an incorrect mass estimate of the subdwarf, for 
example an ELM progenitor is assumed rather than a normal sub¬ 
dwarf. Most objects (60-70 per cent) are low mass main sequence 
stars, while 30-40 per cent are white dwarfs. 

Our own survey delivered a 37 per cent fraction of hot sub¬ 
dwarfs in close binaries (Section 3.1.3). Our survey strategy was 
aimed at and successfully uncovered short-period binaries. Fig.|19| 
shows that setting our detection threshold at 10 km s“* would have 
allowed for the detection of any stellar companion with an orbital 
period < 20 d, any late-type companion with a mass of 0.3 Mq and 
P < 600 d, or just about any white dwarf companion. However, 
a close examination of data sampling shows that 60 per cent were 
obtained with a span of 2 days or less with another 30 per cent 
with a span of 100-400 days, i.e., during a subsequent observing 
run or season. Systems with periods of 10-20 days or longer would 
have only been partially covered and most likely avoided detection. 
Note that the longest period detected in our survey is only 3 d long. 
Setting our detection threshold at 10 d, i.e., some 155 objects out 
of 179 known binaries, the total yield including longer period bina¬ 
ries could be 3! 15 per cent larger for a total binary fraction of 43 per 
cent. The four additional hot subdwarfs with composite spectra (see 
Table [TJ which are likely to have longer orbital periods (> 10 d) are 
such objects. 


4.2 Properties of known binaries: Kinematics 


We calculated the Galactic velocity components (17, V, IV), which 
are relative to the local standard of rest (LSR), of all known hot 
subdwarf binary systems (listed in Appendix C) using their posi¬ 
tions, systemic (y) velocities, proper motions and apparent magni¬ 
tudes. We adopted the right-handed system for the velocity com¬ 
ponents, where U is positive in the direction of the Galactic cen¬ 
tre, V is positive in the direction of Galactic rotation and IV is 
positive toward the North Galactic Pole. We assumed that the so¬ 
lar motion relative to the LSR is ((/, V, IV) = (10.1,4.0,6.7) km s“' 
as determined by |Hogg et al.| ( [^05[ >. The distribution of systemic 
velocities, i.e., radial velocities, follows N ~ where cr^ = 

41 km s“*. The CTr value is the one-dimension equivalent of the two- 
dimensio n transversal velocity dispersion c tt = 59kms“* mea¬ 


sured by Venues, Kawka & Nemeth i2011 1 , where ctj = flcTr. 
In their study of the kinematics of EHB stars, |Altmann et aL]p004| > 
measured significantly larger radial velocities with a distribution 
following CTr = 65kms“* compared to our sample, but they mea¬ 
sured a transversal velocity distribution consistent with the present 


To calculate the Galactic velocity vectors, we employed the 
method outlined in |Johnson & Soderblom|(|1987|l using as a n in¬ 
put the radial velocity, proper motion ( [Zacharias et al.|201^ and 
distance measurements. We determined the distance toward each 
subdwarf using the distance modulus V — My = 5 log d - 5, where 
the magnitude V is listed in Table [CT] We estimated the absolute 
magnitude My from the measured stellar parameters, i.e.. 


My = -2.5 log {AtiQ. Hy), 

where Cl = d with d = 10 pc and d = GMjg. The Eddington 
flux is averaged (Hy) over the Johnson V transmission curve. We 
assumed M = 0.23 Mq for low-mass objects, 0.47 Mq for normal 
objects and 1.0 Mq for subdwarfs in massive binaries, and the pub¬ 
lished surface gravity measurements were usually obtained using a 
spectroscopic method similar to that described in Section 3.1. 

Fig. 1^ shows the U and V velocity components for all hot 
subdwarf binary systems. Table|^lists the ((/, V, IV) velocity com¬ 
ponents and dispersions for the sample of known binaries. The dis¬ 
tribution appears asymmetric with several objects trailing at large 
negative Galactic V velocity. We computed the straight average and 
dispersion (“All”) but excluding the extreme case of OGLE BUL- 
SC16335. We also fitted the distributions with Gaussian functions 
excluding outliers, i.e., all bins with less than three members (‘A > 
3 ‘). The sample velocity dispersion is significantly smaller than that 
calculated by |Altmann et al.|j2004) for single hot subdwarf stars 
{(Ty = 74kms“\ (Ty = 79kms“', CTyy = 64kms“*)- However, 
our velocity dispersion is in better agreement with the dispersion 
{(Ty = 62 ± 8kms“', cry = 52 ± 7kms“', cr^r = 59 ± 8kms“*) 
calculated by |de Boer et al.| ( |1997^ based on a sample of 41 hot 
subdwarf stars. Possible explanations for the inflated Galactic ve¬ 
locities of |Altmann et al.|p004| l are that their sample included yet 
unidentified binaries or that lower dispersion spectroscopy resulted 
in larger measurement errors. In summary, the hot subdwarf pop¬ 
ulation and the confirmed binaries among them are kinematically 
indistinguishable and drawn from the same, general population of 
EHB stars. 

Also, Table [^compares results for the hot subdwarf popula¬ 
tion with that of field white dwarf stars. [Pauli et al.| ( f2006t list 361 
thin disc members and 27 thick disc members, while [Kawka et alTj 
( |2012b[ l list 57 old disc white dwarfs, which is a mix of old thin- 
disk and thick-disk populations, and at least one halo white dwarf. 
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A comparison of the velocity dispersions shows that the binary pop¬ 
ulation may be an even older population than field white dwarf stars 
with some population members belonging to the old disc or even 
the halo. 

Four systems display peculiar kinematics. Three of these, 
SDSS J1622-H4730, PHL 861 and SDSS 11505-^1108, have Galac¬ 
tic velocities that make them halo candidates with a few additional 
objects lagging in the V component making them thick-disk can¬ 
didates. The fourth object, OGLE BUL-SC16335, is in a crowded 
field and there is a possibility that the proper motion measurements 
are incorrect. 

|Barlow et al.| ( |2013a| > calculated kinematics for 5 long pe¬ 
riod systems and found that two of these (PG 1449-1-653 and 
PG 1701-1-359) have kinematics suggesting that they belong ei¬ 
ther to the thick disc or halo. They also report that there is a high 
probability that PG 1104-)-243 belongs to the thick disc. Our cal¬ 
culated Galactic velocities are similar to those of IBarlow et al.l 
) |2013a| >. Note that |Barlow et al.|j2013a| ) includes the disc rotation 
(220 km s“*) in their V velocity. 

A comparison of the velocity components of the hot subdwarf 
binary population to that of |Kawka et al.|j2012b| ) shows that the dis¬ 
persion is larger for all velocity components than that of the white 
dwarf population, suggesting that the subdwarf population appears 
to be older than the white dwarf population. Finally, a comparison 
with the work of [Pauli et al.| j2006) shows that the hot subdwarf 
binary population has a velocity dispersion between the thin and 
thick disc dispersions for white dwarfs. 


4.3 Low mass subdwarfs as progenitors of extremely 
low-mass white dwarfs 


The first known ELM white dwarf progenitor, FID 188112, was dis¬ 
covered by [Heber et al.j (|2003|(. As part of their survey of ELM 
white dwarfe, [Kilic et al.| ) |201lt found that SDSS J1625-H3632 
is similar to FID 188112. Other recently discovered systems are 


KIC 6614501 ( jSilvotti et al.|2012t and NGC 6121-V46 ( [Q’Toolej 
jet al.|[2006) . Our radial velocity survey adds one more object to 
the small sample of ELM white dwarf progenitors :|Vennes, Kawka| 
|& NemeThI j201 1| > showed that GALEX J0805-1058 has atmo¬ 


spheric properties representative of ELM white dwarf progeni¬ 
tors (M <0.3 Mq) and, therefore, it was selected for radial ve¬ 
locity follow-up measurements. New radial velocity measurements 
proved that GALEX J0805-1058 is in a close binary, and, because 
it also lies below the ZAEFIB (see Section 3.1), we conclude that it 
is a genuine ELM white dwarf progenitor. It is also the first ELM 
white dwarf progenitor without a more massive white dwarf com¬ 
panion, and this is likely to have significant implications for the 
origin and evolution of ELM white dwarfs. Two other objects from 
our sample lie below the ZAEHB, J1411-1-7037, which is paired 
with an F star, and J2153-7004, although we did not detect signif¬ 
icant radial velocity variations. 


4.4 Summary 

We presented an analysis of the orbital properties of seven 
new systems comprising a hot subdwarf primary. The secondary 
in one system is a late-type star showing a reflection effect 
(GALEX J2205-3141), while we found evidence that the sec¬ 
ondary star in GALEX J0805-1058 is a very low mass M dwarf 
or possibly a substellar object. The mass function of the other ob¬ 
jects implies the likely presence of a white dwarf companion. The 


Table 4. Kinematical properties of the hot subdwarf binary population. 


Vel." 

N>3'’ 

AID 

sdB'' 

WD^ 

WD,Wn 

WD,h,ck f 

U 

0±5 

2 

-8 

-7.8 



o-u 

52 ±3 

62 

74 

42.8 

34 

79 

V 

-30 ±2 

-32 

-37 

-40.1 


-52 

o-v 

42 ±2 

47 

79 

31.9 

24 

36 

W 

-5 ±3 

-6 

12 

-5.9 



cr w 

34 ±2 

41 

64 

27.4 

18 

46 


All velocities expressed in km s * 
^ This work, but excluding outliers 
This work 


Altmann et 

al. 

2004 

Kawka et al 

i 

012 b 

Pauli et al. 

2006 



- ' ' ' ' 1 ' ' ' ' 1 ' ^ 

''' 1' 

- > 

J 

r 1 

u, 

1 



^ 15 



V (km s->) N 


Figure 21. Galactic velocity vectors U and V of all known binaries con¬ 
taining a hot subdwarf. The individual distributions are shown in the upper 
panel (V) and right panel (17). Details of the measurements are shown in 
Table|4]and discussed in Section 4.2. 


period of photometric variability of two additional systems is, in the 
case of GALEX J1736-1-2806 probably coincident with the orbital 
period, and, in the case of J2038-2657, probably coincident with 
the rotation period of the giant companion. Our survey results, tak¬ 
ing into account the survey strategy, imply an incidence of binarity 
of -“43 per cent in the hot subdwarf population. 

We have compiled a list of all known hot subdwarfs in binary 
systems and performed a binary population analysis. We found that 
systems showing the reflection effect have components that are of 
a lower mass (0.08 to 0.30 Mq) than those that do not show the re¬ 
flection (~0.6Mo). It is very likely that the companion to the hot 
subdwarf in most of the systems not showing the reflection effect in 
the short period binaries (P < 1 d) are white dwarfs. The inferred 
secondary mass distribution is a superposition of two approximate 
power laws, one low-mass power-law (>0.1 Mq) and composed of 
low-mass main-sequence stars, and another, high-mass power-law 
(>0.5 Mo) and primarily composed of white dwarfs with a few 
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early-type main-sequence stars. White dwarfs constitute x!30-40 
per cent of all binary companions. 

We have calculated the Galactic velocity components for all 
known hot subdwarfs in binary system and showed that this pop¬ 
ulation may be older than the field white dwarf population. In this 
sample, we found three systems that possibly belong to the halo. 

Future work will involve high-dispersion spectroscopic 
follow-up of low-velocity amplitude binary candidates, and of bi¬ 
naries comprising a hot subdwarf and an early-type main-sequence, 
or giant companion. 
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APPENDIX A: PHOTOMETRY AND SPECTRAL 
ENERGY DISTRIBUTIONS. 

Fig .|Al[ Fig. |A2| and Fig. |A3| show available photometry from Ta- 
ble |Al [ compared to model spectra. The extinction was determined 
by comparing the observed photometry with the model spectrum. 
Extinction was measured toward several sdB stars (e.g., |Moehler etj 
|al.|1990[|Aznar Cuadrado & Jeffery|2001| l and given that these stars 
are spread through out the Galaxy, i.e., at high latitudes and in the 
Galactic plane, the extinction coefficients E{B — V) can vary from 
0 to as much as 0.4. In these cases the extinction was measured 
using International Ultraviolet Explorer {WE) spectra combined 
with optical and infrared photometry. We used the parameterized 
extinction law as defined by jCardelli et al.| ( |T989^ with R = 3.2. 

We acquired from the Mikulski Archive for Space Telescopes 
(MAST) the following set of WE spectra: 

J0047-1-0958: swp26276mxlo and lwp07169mxlo, 
J0059-I-1544: swp27142mxlo and lwp07144mxlo, 
J1435-)-0013: swp23176mxlo and lwp03501mxlo, 
J1632-)-0759: swp33790mxlo and Iwp 13481 mxlo, 

J1902-5130: swpl7051mxlo and lwrl3321mxlo, 

J2344-3426: swpl7981mxlo. 

The model distributions include the effect of interstellar ex¬ 
tinction with the line-of-sight extinction coefficient obtained from 
[Schlegel, Einkbeiner, & Davis| ( |199^ . Also, we experimented with 
larger coefficients in an attempt to match SEDs showing possible 
intrinsic absorption. For example, the SED of J1632-t-0759 reveals 
the presence of an infrared flux excess as well as a larger extinction 
that revealed injSchlegel, Finkbeiner, & Davis|([T998)’s maps. 
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Table Al. Photometric measurements. 
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Figure Al. Spectral energy distribution of the observed GALEX sample. 


© 2015 RAS, MNRAS 000,[T|t36l 



























f„(erg cm"* s"‘ Hz" 


26 


Kawka et al 


J0751+0925 



J1041-0730 



10 - 8 * 


10-86 


10-88 


io-«' 


J1435+0013 




X(A) 


Figure A2. Spectral energy distribution of the observed GALEX sample. 
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Figure A3. Spectral energy distribution of the observed GALEX sample. 
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APPENDIX B: RADIAL VELOCITY MEASUREMENTS. 

Table |B 1 [ lists the heliocentric-corrected velocity measurements (u), 
the mid-exposure heliocentric Julian dates (HID) and data sources. 
The instrument conhgurations and error estimates are described in 
Sections 2.1 and 2.2, respectively, and the measurement procedures 
are described in Section 3.1. 

Stars labelled with the suffix “B” are the “secondary” compo¬ 
nents of each systems. The radial velocities are those of that sec¬ 
ondary component. 

The label “SSO” refers to spectra obtained with the Wide Field 
Spectrograph (WiFeS) attached to the 2.3m telescope at Siding 
Spring Observatory; the label “KPNO” refers to spectra obtained 
with the 4-m telescope and R.-C. Spectrograph at Kitt Peak Na¬ 
tional Observatory; the label “NTT” refers to spectra obtained us¬ 
ing EFOSC2 attached to the New Technology Telescope at La Silla; 
the label “OND” refers to spectra obtained using the coude spec¬ 
trograph and the 2-m telescope at Ondfejov Observatory; the label 
“SAAO” refers to spectra obtained at the South African Astronom¬ 
ical Observatory using the 74-inch telescope and the Cassegrain 
spectrograph; and new and archival spectra are labelled “FEROS” 
for the Fiber-fed Extended Range Optical Spectrograph on the 
MPG 2.2-m telescope at La Silla, and “INT” and “WHT” for spec¬ 
trographs attached to the Isaac Newton Telescope and William Her- 
schel Telescope at La Palma. 


APPENDIX C: HOT SUBDWARE BINARY SYSTEMS 

We have compiled a list of all known hot subdwarf binary systems. 
Table |Cl| lists the name of the system, its coordinates, proper mo¬ 
tion, V magnitudes, orbital period in days, the systemic velocity (y) 
and the velocity semiamplitude of the hot subdwarf (K). The table 
also lists the eccentricity (e) of the system if the orbit of the binary 
was found to be eccentric. If the companion to the hot subdwarf is 
known it is listed and also if the system is photometrically variable, 
this is also noted in the table. The proper motions and V magnitudes 
for the objects are from The fourth U.S. Naval Observatory CCD 
Astrograph Catalog (UCAC4 [Zacharias et al.|2013) l. Eor those ob¬ 
jects that are not in the UCAC4 and has SDSS photometry, we cal¬ 
culated V magnitudes using the transformation equation of|Jester| 

|eraL| ( |2005t : 

V = g - 0.58(g - r) - 0.01 

The references for the binary system properties are provided in 
the final column. GALEX J1736+2806 and J2038-2657 are not 
included in Table [CT] because the photometric variations are not 
clearly associated to orbital periodicities. New Kepler identifica¬ 
tions (KIC) are added to the list despite current lack of radial veloc¬ 
ity data because of the higher quality of light curve analysis and the 
timeliness of the results. Table |C2| hsts published stellar parameters 
(with relevant references) and calculated absolute V magnitude and 
Galactic velocity vectors (U, V, W). The mass of sdO companions 
to Be stars is assumed to be 1 Mq and that of low-mass sdB stars is 
assumed to be 0.23 Mq. 
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Table Bl. Radial velocities. 


HJD 

(2450000+) 

V 

(kms-') 

Source 

HJD 

(2450000+) 

V 

(kms^*) 

Source 

HJD 

(2450000+) 

V 

(kms^*) 

Source 

HJD V 

(2450000+) (kms-‘) 

Source 

GAT.F.X J0047+0337B 

A = 4, D = 

: 45.3, Cr„ : 

= 4.4 

5930.74878 

31.3 

KPNO 

6047.05680 29.3 

sso 

5898.04458 

-19.9 

sso 

GALEX J0401-3223 

5930.80712 

18.8 

KPNO 

6047.08829 31.4 

sso 

5899.04507 

-23.4 

sso 

3953.91980 

53.0 

EEROS 

5930.86267 

13.6 

KPNO 

6700.41812 22.9 

SAAO 

5930.57007 

-22.1 

KPNO 

4014.75477 

52.1 

EEROS 

5930.91715 

5.7 

KPNO 

6700.47565 34.4 

SAAO 

5930.57406 

-23.5 

KPNO 

5500.79501 

49.3 

EEROS 

5931.00340 

0.5 

KPNO 

6700.52813 30.5 

SAAO 

5930.66341 

-23.2 

KPNO 

5900.02546 

45.2 

SSO 

5932.71032 

50.9 

KPNO 

A = 6, S = 29.8, cr„ = 

3.5 

5931.56622 

-26.2 

KPNO 

6171.78071 

55.8 

NTT 

5932.86851 

42.3 

KPNO 

GALEX J1114-2421 

5932.55868 

-27.1 

KPNO 

6171.83654 

66.3 

NTT 

5932.95043 

36.1 

KPNO 

6046.17066 20.2 

SSO 

6172.84455 

-5.7 

NTT 

6172.75132 

56.8 

NTT 

5933.00770 

30.1 

KPNO 

6047.00935 19.7 

sso 

M 

oo' 

II 

: -21.4, (T„ 

= 6.3 

6172.78790 

52.3 

NTT 

A = 9, S = 

25.5, (T„ = 

16.0 

6047.04136 15.7 

sso 

GALEX J0047+0958 (HD 4539) 

6173.74065 

72.0 

NTT 

GALEX J0751+0925 

6047.10305 15.9 

sso 

3227.46321 

3.1 

OND 

6173.80204 

63.0 

NTT 

5898.15513 

164.5 

SSO 

6047.99465 16.9 

sso 

3227.48556 

-0.4 

OND 

6173.88283 

64.4 

NTT 

5898.20329 

-13.8 

SSO 

6700.43054 7.6 

SAAO 

3227.50802 

-1.4 

OND 

6700.28962 

54.0 

SAAO 

5899.12739 

-113.9 

SSO 

6700.48786 4.8 

SAAO 

3229.51947 

-2.4 

OND 

6700.34266 

49.4 

SAAO 

5899.22915 

170.8 

sso 

6700.53977 16.3 

SAAO 

3229.54285 

-0.2 

OND 

6700.40330 

46.9 

SAAO 

5930.75957 

87.8 

KPNO 

= 8, V = 14.6, (T„ = 

5.2 

3255.46637 

-0.5 

OND 

A = 14, S = 

= 55.8, cr„ 

= 7.6 

5930.81729 

92.2 

KPNO 

Feige 66 


3255.49221 

-2.4 

OND 

GALEX J0500+0912 

5930.87413 

-127.6 

KPNO 

5311.46603 -5.5 

OND 

3255.54837 

-3.0 

OND 

5930.62297 

48.5 

KPNO 

5930.92836 

42.6 

KPNO 

5311.48899 -7.9 

OND 

3255.57093 

-2.1 

OND 

5930.73469 

47.7 

KPNO 

5930.93931 

96.5 

KPNO 

5311.51191 -4.2 

OND 

5499.73722 

-3.2 

EEROS 

5931.65783 

50.2 

KPNO 

5930.98374 

135.2 

KPNO 

5311.53487 5.8 

OND 

5500.56785 

-3.7 

EEROS 

5932.74346 

45.0 

KPNO 

5931.03215 

-91.9 

KPNO 

5311.55792 -7.1 

OND 

5500.70433 

-2.7 

EEROS 

A = 4, D = 

: 47.8, 0-„ : 

= 1.9 

5931.73038 

-8.1 

KPNO 

5312.46039 -3.2 

OND 

5500.78859 

-3.2 

EEROS 

GALEX J0507+0348 

5931.87541 

141.4 

KPNO 

5312.48335 -5.0 

OND 

5930.55564 

-5.0 

KPNO 

5931.67914 

77.7 

KPNO 

5931.99730 

41.9 

KPNO 

5312.50638 0.0 

OND 

5930.55774 

0.0 

KPNO 

5931.83227 

39.7 

KPNO 

5932.75723 

163.1 

KPNO 

5312.52952 -4.2 

OND 

5930.65685 

-3.4 

KPNO 

5931.92094 

91.2 

KPNO 

5932.88029 

-4.4 

KPNO 

5312.55266 -4.7 

OND 

5931.57524 

-3.7 

KPNO 

5932.69831 

116.5 

KPNO 

5933.01902 

-125.4 

KPNO 

= 10, V = -3.6, = 

= 3.7 

5932.55265 

-6.4 

KPNO 

5932.83768 

30.7 

KPNO 

6046.91015 

60.7 

SSO 

GALEX J1356-4934 

6172.83793 

0.1 

NTT 

5932.89420 

28.6 

KPNO 

6047.93820 

163.7 

SSO 

6046.15529 0.2 

SSO 

6173.77623 

-1.4 

NTT 

6171.78863 

147.5 

NTT 

A = 19, S = 

46.1, (T„ = 

101.1 

6046.27506 17.2 

sso 

N = 20, V 

= -2.1, cr„ 

= 2.1 

6171.84440 

155.4 

NTT 

GAT.F.X .10805-1058 

6047.13503 4.8 

sso 

GALEX J0049+2056 

6171.85811 

153.7 

NTT 

5898.12316 

52.8 

SSO 

6172.47898 12.8 

NTT 

5930.59893 

15.1 

KPNO 

6171.90134 

154.1 

NTT 

5898.12983 

64.7 

sso 

6172.49297 7.1 

NTT 

5930.68516 

17.7 

KPNO 

6172.81341 

140.7 

NTT 

5898.13875 

69.2 

sso 

6172.51996 15.8 

NTT 

5931.58588 

14.6 

KPNO 

6172.90675 

166.4 

NTT 

5898.25067 

26.9 

sso 

6172.54716 8.9 

NTT 

5932.57896 

18.2 

KPNO 

6173.81997 

99.9 

NTT 

5899.08066 

61.3 

sso 

6700.49865 20.6 

SAAO 

N = 4, v: 

= 16.4, cr„ : 

= 1.6 

6173.90692 

166.6 

NTT 

5899.20870 

88.4 

sso 

6700.56149 8.7 

SAAO 

GALEX J0059+1544 

6700.31140 

58.5 

SAAO 

5900.24103 

78.0 

sso 

A = 9, S = 10.7, cr„ = 

6.1 

3955.82658 

20.6 

EEROS 

6700.36337 

111.9 

SAAO 

5930.76915 

47.5 

KPNO 

GALEX J1407+3103 

3986.81551 

9.6 

EEROS 

A= 16, D = 

108.7, o-„ 

= 47.9 

5930.82802 

88.9 

KPNO 

5932.03194 -160.6 

KPNO 

5930.58194 

15.5 

KPNO 

GALEX J0613+3420 

5930.88403 

42.1 

KPNO 

5932.97687 -157.3 

KPNO 

5930.58594 

17.3 

KPNO 

5931.74486 

-135.5 

KPNO 

5930.97349 

75.2 

KPNO 

5933.05659 -164.7 

KPNO 

5930.67212 

20.7 

KPNO 

5931.84915 

-139.9 

KPNO 

5930.99445 

92.6 

KPNO 

A = 3, S = -160.9, o-„ 

= 3.0 

5931.60011 

13.1 

KPNO 

5931.93741 

-138.7 

KPNO 

5931.75141 

42.1 

KPNO 

GALEX J1411+7037B 

5932.56969 

14.9 

KPNO 

5932.72107 

-40.9 

KPNO 

5931.88363 

85.1 

KPNO 

5932.00693 -18.0 

KPNO 

5932.72419 

18.3 

KPNO 

A = 4, D = - 

•113.7, cr„ 

= 42.1 

5931.98199 

45.6 

KPNO 

5932.04748 -17.6 

KPNO 

6172.83077 

12.4 

NTT 

GALEX J0657-7324 

5932.86147 

49.1 

KPNO 

5932.99447 -14.9 

KPNO 

6173.78258 

21.6 

NTT 

6047.87067 

-7.0 

SSO 

5932.94370 

73.5 

KPNO 

5933.06557 -20.9 

KPNO 

N=l0,v 

= 16.4, cr„ 

= 3.8 

6171.82808 

-11.3 

NTT 

6046.92589 

31.8 

SSO 

A = 4, ii = -17.9, cr„ = 

= 2.1 

GALEX J0206+1438 

6171.92035 

-17.6 

NTT 

6046.96533 

36.0 

SSO 

GALEX J1421+7124 

5932.64012 

19.8 

KPNO 

6172.80377 

-11.7 

NTT 

6047.95401 

48.4 

SSO 

5461.35716 -24.2 

OND 

5932.76929 

23.2 

KPNO 

6172.91817 

1.1 

NTT 

6700.39670 

36.6 

SAAO 

5461.38072 -22.3 

OND 

6172.76462 

19.9 

NTT 

6172.92477 

-12.6 

NTT 

6700.46184 

59.2 

SAAO 

5461.40360 -25.5 

OND 

6172.79845 

4.1 

NTT 

6700.33004 

-16.0 

SAAO 

6700.51455 

86.2 

SAAO 

5462.34435 -29.5 

OND 

6172.89971 

9.8 

NTT 

6700.37998 

-13.9 

SAAO 

A = 23, D = 

: 60.1, (T„ = 

20.0 

5462.36727 -19.9 

OND 

6173.75136 

1.6 

NTT 

6700.44564 

-19.7 

SAAO 

GALEX J0812+1601 

5462.64138 -17.4 

OND 

6173.81315 

12.5 

NTT 

II 

SO 

Cl 

II 

-12.1, o-„ 

= 5.8 

5898.18568 

-26.7 

SSO 

5463.32822 -21.2 

OND 

6173.92187 

19.6 

NTT 

GALEX J0703+6236 

5898.23480 

-20.9 

sso 

5463.35791 -25.7 

OND 

N = S, v : 

= 13.8, o-„ : 

= 7.5 

5931.04474 

24.2 

KPNO 

5932.85141 

1.3 

KPNO 

5464.32919 -23.3 

OND 

GALEX J0232+4411 

5931.71772 

18.3 

KPNO 

5932.93361 

-1.6 

KPNO 

5464.35893 -24.2 

OND 

5931.61405 

41.1 

KPNO 

5931.85802 

17.6 

KPNO 

5933.02934 

-4.1 

KPNO 

5470.25157 -24.4 

OND 

5931.81305 

52.1 

KPNO 

5931.94502 

19.8 

KPNO 

A = 5, « = 

-10.4, cr„ = 

: 11.2 

5483.25104 -35.4 

OND 

5932.65498 

46.2 

KPNO 

A = 4, D = 

: 20.0, (T„ : 

= 2.6 

GAT,EX J1041-0730 

5794.51260 -28.6 

OND 

5932.81723 

41.8 

KPNO 

GALEX J0716+2319 

6046.11527 

30.0 

SSO 

A= 13, D = -24.7, o-„ 

= 4.4 
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Table Bl. continued 


HJD 

(2450000+) 

V 

(kms-') 

Source 

HJD 

(2450000+) 

V 

(km s“*) 

Source 

HJD 

(2450000+) 

V 

(kms^*) 

Source 

HJD 

(2450000+) 

V 

(kms^*) 

Source 

GALEXJI427-2701 

5757.01406 

-4.4 

sso 

5757.09665 

7.1 

sso 

5898.96627 

66.8 

sso 

6048.11285 

-6.6 

sso 

6048.21613 

52.3 

SSO 

5757.10547 

14.4 

sso 

5898.99478 

67.3 

sso 

6171.47061 

10.9 

NTT 

6171.57012 

28.2 

NTT 

5760.12396 

14.0 

sso 

5899.02602 

60.4 

sso 

6171.48306 

6.8 

NTT 

6171.59018 

27.1 

NTT 

5761.03781 

20.1 

sso 

5899.98304 

48.7 

sso 

6171.49578 

1.1 

NTT 

6171.60924 

17.8 

NTT 

5761.09468 

8.6 

sso 

6171.72238 

-38.2 

NTT 

6171.50953 

3.5 

NTT 

6171.63315 

13.7 

NTT 

6171.68438 

10.4 

NTT 

6171.77618 

-50.8 

NTT 

6171.52366 

5.2 

NTT 

6172.56912 

51.2 

NTT 

6171.74097 

10.7 

NTT 

6171.82488 

-69.0 

NTT 

6171.53666 

-0.5 

NTT 

6172.60585 

42.9 

NTT 

6172.65633 

18.1 

NTT 

6171.89529 

-78.4 

NTT 

6172.50574 

0.6 

NTT 

6172.62568 

35.3 

NTT 

6172.72634 

4.9 

NTT 

6172.63823 

72.5 

NTT 

6172.53333 

-7.6 

NTT 

6361.86864 

-70.0 

NTT 

N = 9,v-- 

= 12.0, cr„ = 

= 4.7 

6172.70799 

50.3 

NTT 

6700.57282 

-12.0 

SAAO 

6361.90960 

-84.9 

NTT 

GAT .EX J1903-3528 

6172.77525 

28.4 

NTT 

6700.60630 

-5.0 

SAAO 

6429.53360 

51.7 

INT 

6171.70234 

-29.7 

NTT 

6172.85901 

-11.3 

NTT 

A= 11, D 

II 

1 

o 

= 6.6 

6870.49571 

-17.8 

NTT 

6171.75097 

-37.6 

NTT 

6173.67386 

71.5 

NTT 

GALEXJ1435+0013 

6870.61105 

29.4 

NTT 

6172.67409 

-17.3 

NTT 

6173.73606 

98.9 

NTT 

3829.84393 

-5.5 

EEROS 

6870.69156 

29.8 

NTT 

6172.73670 

-25.6 

NTT 

6173.79014 

96.8 

NTT 

4137.79323 

-3.2 

EEROS 

N=16,v = 

11.4, o-„ = 

= 40.0 

6173.70238 

-2.9 

NTT 

6173.87812 

55.8 

NTT 

5756.90232 

-16.0 

SSO 

GALEXJ1736+2806B 

N = 5,d = 

-22.6, (T„ : 

= 11.8 

6870.63480 

96.2 

NTT 

6046.18707 

-3.7 

SSO 

5045.37797 

-14.2 

OND 

GAT.EX .11911-1406 

6870.71233 

93.8 

NTT 

6047.22908 

-6.5 

SSO 

5045.39926 

-16.5 

OND 

5760.14357 

-156.8 

SSO 

6870.86537 

58.6 

NTT 

6047.26718 

2.2 

SSO 

5310.44503 

-19.6 

OND 

5761.05696 

-147.0 

sso 

6989.54548 

54.3 

EEROS 

6048.04105 

3.6 

SSO 

5310.46808 

-8.2 

OND 

5761.12734 

-151.3 

sso 

6990.55028 

-34.2 

EEROS 

6048.07142 

-5.9 

SSO 

5310.49159 

-16.9 

OND 

N = 3,v = 

-151.7, o-„ 

= 4.0 

6994.58156 

75.6 

EEROS 

6048.09997 

-4.6 

SSO 

5310.51449 

-7.4 

OND 

GAT .EX .12153-7003 

N = 24, V-- 

= 32.6, (T„ : 

= 55.0 

6700.58473 

-12.6 

SAAO 

5310.53754 

-14.7 

OND 

3918.75562 

43.4 

EEROS 

GALEX J2334+5347 

6700.62607 

-6.0 

SAAO 

5310.56068 

-9.3 

OND 

3956.76278 

42.5 

EEROS 

5461.43461 

40.3 

OND 

A= 11, D 

II 

1 

= 5.3 

5310.58391 

-13.9 

OND 

5757.17717 

41.8 

SSO 

5461.48740 

36.1 

OND 

J1600-6433 


5311.58102 

-13.2 

OND 

5757.18600 

43.1 

SSO 

5461.55799 

31.5 

OND 

5756.93811 

55.6 

SSO 

5311.60395 

-18.7 

OND 

5761.21680 

44.2 

SSO 

5461.58084 

44.7 

OND 

5756.94716 

55.1 

sso 

5312.57540 

-11.5 

OND 

5897.97485 

41.4 

SSO 

5462.39976 

46.5 

OND 

5757.03853 

48.6 

SSO 

5312.59830 

-7.9 

OND 

5897.98518 

36.9 

SSO 

5462.42295 

41.0 

OND 

5757.04736 

53.6 

sso 

5377.40465 

-11.8 

OND 

5898.94745 

39.8 

SSO 

5462.59518 

31.6 

OND 

5757.96919 

65.3 

sso 

5377.42598 

-6.9 

OND 

5898.95287 

42.1 

SSO 

5462.61797 

42.5 

OND 

5757.97801 

56.6 

sso 

5430.35500 

-9.3 

OND 

5898.97969 

40.3 

SSO 

5463.41196 

37.9 

OND 

5760.98197 

41.6 

sso 

5430.37835 

-10.7 

OND 

5899.96199 

41.4 

SSO 

5463.44187 

29.6 

OND 

6171.62174 

37.8 

NTT 

5483.27549 

-7.8 

OND 

6171.80243 

50.5 

NTT 

5463.47292 

26.0 

OND 

6171.66370 

45.9 

NTT 

5675.52280 

-16.0 

OND 

6171.87909 

47.0 

NTT 

5464.41610 

37.2 

OND 

6172.58713 

56.0 

NTT 

5766.38926 

-7.9 

OND 

6172.68364 

42.9 

NTT 

5464.44582 

50.4 

OND 

6172.66055 

40.6 

NTT 

5794.39180 

-6.5 

OND 

6173.79603 

54.2 

NTT 

5483.37624 

28.3 

OND 

6700.59333 

43.1 

SAAO 

N = 2l,d = 

-11.9, o-„ 

= 4.0 

N=15,v 

= 43.4, (T„ 

= 4.2 

5483.56518 

29.9 

OND 

6700.63352 

37.7 

SAAO 

GALEXJ1753-5007B 

GALEX J2205-3141 

5794.58849 

34.8 

OND 

N=13, V 

= 49.0, o-„ 

= 8.3 

5757.06955 

-56.9 

SSO 

5757.23298 

-67.4 

SSO 

5834.29578 

29.0 

OND 

GALEXJ1632+0759 

5757.07837 

-56.4 

SSO 

6171.77057 

-17.9 

NTT 

A = 17, « 

= 36.3, (T„ 

= 6.9 

5435.41483 

-67.4 

WHT 

5761.02112 

-65.7 

SSO 

6171.81889 

-68.1 

NTT 

GALEX J2344-3426 

5756.97348 

-62.5 

SSO 

5761.07352 

-60.3 

sso 

6171.88946 

-62.9 

NTT 

3679.56424 

21.5 

EEROS 

5761.00101 

-50.0 

SSO 

6171.64403 

-56.2 

NTT 

6172.70173 

23.7 

NTT 

3957.88130 

21.5 

EEROS 

6048.12788 

15.2 

SSO 

6171.71118 

-73.0 

NTT 

6172.78429 

-8.0 

NTT 

5757.29066 

12.9 

SSO 

6171.55953 

-0.7 

NTT 

6172.61645 

-75.6 

NTT 

6172.85277 

-62.3 

NTT 

5898.02545 

21.5 

SSO 

6171.57889 

8.6 

NTT 

6173.71129 

-74.9 

NTT 

6173.73046 

32.9 

NTT 

5899.01277 

19.2 

SSO 

6172.55735 

4.8 

NTT 

A = 8, u = 

-64.9, cr„ 

= 8.0 

6173.76269 

19.9 

NTT 

5899.06162 

21.0 

SSO 

6172.57774 

-3.1 

NTT 

GALEXJ1845-4138 

6173.83761 

-38.8 

NTT 

6171.81182 

29.0 

NTT 

6432.47442 

-31.5 

INT 

6171.65588 

-52.3 

NTT 

6870.62358 

-49.9 

NTT 

6171.91775 

26.4 

NTT 

6870.48215 

-68.0 

NTT 

6171.73152 

-52.3 

NTT 

6870.70432 

-42.0 

NTT 

6172.69358 

31.2 

NTT 

6870.59970 

-69.0 

NTT 

6172.64693 

-57.2 

NTT 

6870.85814 

14.7 

NTT 

6173.76988 

34.9 

NTT 

6870.68037 

-65.0 

NTT 

6172.71683 

-69.1 

NTT 

A= 13, S = 

: -25.1, (T„ 

= 36.4 

A = 10, S 

= 23.9, o-„ 

= 6.1 

II 

(N 

II 

: -32.4, o-„ 

= 33.3 

6173.68314 

-56.9 

NTT 

GALEX J2254-5515 




GALEXJ 1731+0647 

N = 5.v = 

-57.6, cr„ 

= 6.1 

5898.00341 

-7.3 

SSO 




5757.00524 

-19.9 

SSO 

GALEXJ 1902-5130 

5898.05579 

-24.7 

sso 
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Table Cl. Properties of hot subdwarf binary systems. 
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Table C2. Kinematics (U, V, IK), stellar parameters and absolute V magnitude (My) of known binaries. 
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